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Abstract. We present the results of analysis of "snapshot" spectra of 253 metal-poor halo stars —3.8 < [Fe/H] < —1.5 
obtained in the HERES survey. The snapshot spectra have been obtained with VLT/UVES and have typically S/N ~ 54 per 
pixel (ranging from 17 to 308), R ~ 20000, A = 3760-4980 A. This sample represents the major part of the complete HERES 
sample of 373 stars; however, the CH strong content of the sample is not dealt with here. 
The spectra are analysed using an automated line profile analysis method based on the Spectroscopy Made Easy (SME) codes 
of Valenti & Piskunov. Elemental abundances of moderate precision (absolute rms errors of order 0.25 dex, relative rms errors 
of order 0.15 dex) have been obtained for 22 elements, C, Mg, Al, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Zn, Sr, Y, Zr, Ba, La, 
Ce, Nd, Sm, and Eu, where detectable. Of these elements, 14 are usually detectable at the 3a confidence level for our typical 
spectra. The remainder can be detected in the least metal-poor stars of the sample, spectra with higher than average S/N, or 
when the abundance is enhanced. 
Among the sample of 253 stars, disregarding four previously known comparison stars, we find 8 r-II stars and 35 r-I stars. 
The r-II stars, including the two previously known examples CS 22892-052 and CS 31082-001, are centred on a metallicity 
of [Fe/H] = —2.81, with a very small scatter, on the order of 0.16 dex. The r-I stars are found across practically the entire 
metallicity range of our sample. We also find three stars with strong enhancements of Eu which are s-process rich. A significant 
number of new very metal-poor stars are confirmed: 49 stars with [Fe/H] < —3 and 181 stars with —3 < [Fe/H] < —2. We 
find one star with [Fe/H] < —3.5. 

We find the scatter in the abundance ratios of Mg, Ca, Sc, Ti, Cr, Fe, Co, and Ni, with respect to Fe and Mg, to be similar 
to the estimated relative errors and thus the cosmic scatter to be small, perhaps even non-existent. The elements C, Sr, Y, 
Ba and Eu, and perhaps Zr, show scatter at [Fe/H] £ —2.5 significantly larger than can be explained from the errors in the 
analysis, implying scatter which is cosmic in origin. Significant scatter is observed in abundance ratios between light and heavy 
neutron-capture elements at low metallicity and low levels of r-process enrichment. 
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1. Introduction 

The Hamburg/ESO R-process Enhanced Star (HERES) sur- 
vey has been described in the preceding paper in this series 
(Christlieb et al. 120041 hereafter Paper I). HERES is mostly 
based on confirmed metal-poor stars from the Hamburg/ESO 
survey (HES; Wisotzki et al. l2000l In HERES, "snapshot" 
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spectra of 373 very metal-poor stars, here meaning with 
[Fc/H] < —1.5 1 as judged from medium resolution spectra, 
have been obtained with VLT2-UVES, with the main goal of 
finding stars enhanced in the r-process elements through de- 
tection of strong Eu II lines. Though the snapshot spectra are 
of what would generally be considered low quality for abun- 
dance analysis (typically S/N ~ 54, R ~ 20000, A = 3760- 
4980 A), they contain a wealth of information and abundances 
may be derived for a significant number of elements with mod- 
erate precision (absolute rms errors of order 0.25 dex, rela- 
tive errors of order 0.15 dex). Modern surveys of metal-poor 
stars, such as the ESO "First Stars" Large programme (Cayrel 
et al. 120041 Hill et al. |2002 and refer ences therein), now ob- 
tain significantly better quality spectra for of the order of 70 
stars, yet just a decade ago spectra of similar quality to our 
snapshot spectra were typical for studies of very metal-poor 
stars (e.g. McWilliam et al. I1995al I1995b> . While high pre- 
cision (better than 0.1 dex) is expected to be needed to dis- 
cern detailed patterns in abundance distributions which might 
serve as diagnostics of early nucleosynthesis (e.g. Karlsson & 
Gustafsson 2001), the large number of stars observed in the 
HERES survey offers the possibility to investigate more gen- 
eral trends in metal-poor star abundances in a previously unex- 
plored statistical regime. In particular, the scatter in abundance 
distributions may provide insights into mixing and the diversity 
of supernovae at early epochs. The study of Norris et al. ( 2001 1, 
which investigated such scatter, drawing from different sur- 
veys in the literature, had of order 70 stars in this metallicity 
regime. This project provides a homogeneously analysed sam- 
ple of several hundred stars. 

In this paper we analyse a total of 253 of the spec- 
tra using an automated spectrum analysis technique based 
on the Spectroscopy Made Easy (SME) codes by Valenti & 
Piskunov J1996i . In Sect. |2] the sample, observations, and 
choice of initial stellar parameters will be described. In Sect. [3] 
the automated analysis technique and codes are described, in- 
cluding the line list, detection classification and error analysis, 
and the method is tested through comparisons with previously 
studied stars and tests of the robustness of the method. The 
results are presented in Sect. |4] including the new interesting 
stars found in the survey. In Sect. [5] we discuss general trends 
in the results and other interesting features to emerge. Finally, 
in Sect.[6]we present our conclusions and discuss possible uses 
for the data set. 

2. Observations and Sample 

The target selection for HERES and observational details have 
been described in Sect. 3 of Paper I. The spectra were obtained 
with the ESO-VLT2 and UVES, and cover a wavelength range 
of 3760-4980 A, and have an average signal-to-noise ratio of 
S/N ~ 54 per pixel over the spectral range, though some spec- 
tra have S/N as low as 17 and as high as 308. A 2" slit is 
employed giving a minimum resolving power of R « 20000, 
though typically the resolving power is seeing limited and thus 

1 [ A/B] - log ( N A /N b ) * - log ( N A /N B ) where jV x are number 
densities. 
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Fig. 1. The T c g versus log g Kiel diagram for the sample. The 
diameter of the circle for each star is proportional to [C/Fe]; 
cases where C is undetected are shown with an asterisk. The 
mean error bar is shown in the bottom right. 

slightly better. As mentioned in Paper I, the pipeline-reduced 
spectra are employed, corrected to the stellar rest frame. The 
sample of analysed stars is given in the electronic Tabled 2 with 
coordinates and barycentric radial velocities. The final stellar 
parameters of the sample are presented in Table [21 3 together 
with the abundances for convenience, and are plotted in Fig.^ 
The [Fe/H] distribution for the sample from the final analysis is 
shown in Fig. [2] Though observations of adequate quality have 
been obtained for 373 stars, not all are analysed in this paper. 
Spectra showing strong molecular carbon features cannot be 
analysed by our current method, and thus 72 such spectra from 
the survey are not included in this sample for this reason. These 
spectra are, however, of great interest and it is planned that they 
will be analysed separately in future. Further, some stars which 
turned out to be too Fe-rich ([Fe/H] > —1.5) or too cool for 
our analysis method (< 4200 K), or are suspected to be spec- 
troscopic binaries or rotators, were also removed from the sam- 
ple. For a small number of stars, we have not yet obtained ade- 
quate photometry and they are also removed from the sample. 
Thus, the final sample of stars analysed in this work presented 
in Tables^andlJlcontains 253 stars. 

2.1. Photometry and Effective Temperatures 

To begin our analysis we require as accurate as possible esti- 
mates of the effective temperatures of HERES targets. This is 
particularly the case since, as noted below, the snapshot spectra 
we obtain are not generally of sufficiently high quality to de- 
rive precise spectroscopic temperature estimates. Hence, over 
the course of the past few years, we have obtained broadband 
BVRcIc (where the subscript "C" indicates the Cousins sys- 
tem) observations for as many HERES targets as possible, us- 
ing the ESO/Danish 1 .5m telescope on La Silla and the DFOSC 
instrument. The observing and reduction techniques for this 

2 TableQin its entirety is available only electronically 

3 Table|2|in its entirety is available only electronically 
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Table 1. The sample. The entire table is available only electronically. Coordinates of the HES stars (prefix HE) have been derived 
from the Digitized Sky Survey I and are accurate to 1"; the coordinates of the HK survey stars (prefix CS) are from identifications 
of the sources in the 2MASS All Sky Release. Barycentric radial velocities w rac j were measured from the snapshot spectra and 
are accurate to a few km s _1 . 
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Table 2. Derived stellar atmosphere parameters and elemental abundances for the sample. The entire table is available only 
electronically. This portion of the table is given as a guide to its form. For each star we report the average signal-to-noise ratio 
per pixel for the whole observed spectrum, the stellar parameters with their respective relative and absolute rms error estimates, 
er rcl and er abs , then for each element X the abundance log ex and its relative and absolute rms error estimates (J\ og e x , [X/Fe] 
and its relative and absolute rms error estimates crx/Fel- N gives the number of features of the element used in each star, noting 
that this varies from star to star as certain features may be automatically rejected if the star has a strong G band, if the feature is 
near a Balmer line in warmer stars (see Sect. 13. 2> . or if they are determined to be affected by a cosmic ray hit or bad pixel. N3 
then gives the number of those features classified as 3er detections. In the case of CH bands we simply classify detections and 
non-detections with 1 and respectively. In the case of a non-detection, this is signified by an abundance of —9.99 and error of 
0.00. 
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data are described in Beers et al. (2005, in preparation). We also 
make use of near-infrared JHK photometry from the 2MASS 
Point Source Catalog (Cutri et al. 12003b . The complete set of 
available photometry for stars in the HERES sample will be 
published in Beers et al. (2005, in preparation), a compila- 
tion of photometry for more than 1500 metal-poor stars and 
horizontal-branch candidates. The adopted reddening is taken 
from the maps of Schlegel et al. J1998l >. adjusted to account 



for distance. Note, all but the brightest giants of our sample lie 
outside the reddening layer and thus for most stars the full red- 
dening is applied. The reddening data will also be presented in 
Beers et al. 

Temperature estimates are obtained on the scale of Alonso 
et al. J1999l >. Note that, while the majority of the colour calibra- 
tions of Alonso et al. require measurements on the Johnson sys- 
tem, the near-IR colours J — H and J — K are on the TCS sys- 
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Fig. 2. The distribution of analysed stars in [Fe/H]. The [Fe/H] 
values are from the final analysis of the spectra. 

tem (the photometric system at the 1.54 m Carlos Sanchez tele- 
scope; Arribas & Martinez-Roger [T987> . Thus we first have to 
apply several transformations of our observed colours. We fol- 
low the prescription described by Sivarani et al. (2004i. In per- 
forming our transformations, and for carrying out the Alonso 
et al. estimates of effective temperature, we made use of a pro- 
gram kindly provided to us by T. Sivarani. In order to arrive 
at a final "best estimate" of the effective temperatures for our 
HERES pilot sample stars, we take a straight average of the de- 
rived estimates from different colour criteria, after trimming 
off the highest and lowest estimates for each star. The esti- 
mated effective temperatures, for the colours we had available, 
are provided in Table 13 Although errors in the determinations 
of temperature for individual colours arising from photometric 
errors can range from 50-100 K (for the optical colours), and 
up to several hundred K for the near-IR colours (owing to the 
generally larger errors in the 2MASS photometry), we conser- 
vatively estimate that our final determination of effective tem- 
perature has an absolute rms error of 100 K. We note that our 
photometry is incomplete, i.e. we do not have measurements 
of all colours for every star, and this incompleteness, together 
with the trimming in the averaging procedure, means there is 
an error introduced due to differences in the temperature scales 
from different colour criteria. Since the incompleteness is ef- 
fectively random, this should lead only to a slight increase in 
random scatter. This has been considered in our error estimates. 
The error due to uncertain reddening (see discussion below) 
was also considered. We plan to obtain improved photometric 
estimates of T c g for stars of particular interest on the HERES 
program (for instance, those noted to be r- and/or s-process en- 
hanced) by measuring more precise JHK photometry in the 
near future. 

As will be discussed in Sect. 13.5.11 our temperatures are 
found to be systematically warmer than those found in the lit- 
erature for the stars where comparisons can be made. This dis- 
crepancy was traced primarily to the use of different reddening 
maps. In this work, the maps of Schlegel et al. ( 1998 ) have been 
employed, while in the past most workers have used the maps 




0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 
E(B-V) (Schlegel et al.) 




4000 4500 5000 5500 6000 6500 7000 
T_.„ [Kl 



Fig. 3. Change in E(B — V) (upper panel) and T e g (lower 
panel) if we adopt the reddening maps of Burstein & 
Heiles 



; 1982 1 instead of those of Schlegel et al. (1998). The 
reddening corrections are usually reduced and thus the temper- 
atures become somewhat cooler. 



of Burstein & Heiles ( 1982 1. We computed reddenings and ef- 
fective temperatures using both maps and find systematic dif- 
ferences as seen in Fig. [5] the results using the maps of Burstein 
& Heiles being cooler than those found using Schlegel et al. by 
70 ± 58 K. We have chosen to use the Schlegel et al. estimates, 
which have superior spatial resolution and are thought to have 
a better determined zero point. We note, however, that Arce 
& Goodman ( 1999 1 have pointed out that the Schlegel et al. 
maps may overestimate the reddening values when E(B — V) 
exceeds about 0.15 mag, though none of the stars considered 
here exceed that value (see also Beers et al. 2002 for further 
discussion). 

It should be noted that while the HERES target selection 
aims at a cutoff of B — V > 0.5, a number of warmer stars 
which turned out to be bluer than our targeted cutoff have also 
been included due to incorrect estimated B — V colours from 
the HES prism plates. 

As will be described later in Sect. 13.31 we compute indi- 
vidual abundances for all Fe lines and thus can compare our 
adopted temperature scale with that implied by excitation equi- 
librium. Given the lines available to us, the spectroscopic data 
are not of sufficient quality to determine precise excitation tem- 
peratures, but one can make a comparison of the average tem- 
perature scale. Figure|4]plots the histogram of slopes of the best 
fit to the Fe I line abundances against excitation. On average 
the distribution is skewed towards negative slopes, indicating 
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Fig. 4. Histogram of the slopes of fits to the trends of Fe I line 
abundances with excitation for the sample. 

that our temperature scale is on average warmer than would be 
obtained from excitation equilibrium. Such trends have been 
noted in previous works on metal-poor giants and dwarfs, for 
example by Norris et al. J200U and Cohen et al. J2004> . 
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2.2. Other Atmospheric Parameters 

Estimates of the metallicity, [Fe/H], which serve as initial 
guesses for our automated analysis, are derived from the 
medium-resolution spectra obtained during the course of the 
large-scale campaign to identify suitable targets for the HERES 
program (see Paper I), in combination with available B— V and 
J — K photometry. Estimates of [Fe/H] were obtained from 
the calibration of the Ca II K-line index, KP, along with the 
B — V colour, described by Beers et al. ( 1999 ). For most stars 
this calibration should be accurate to on the order of 0.2 dex, 
given spectra and photometry of suitable quality. In addition, 
we carried out a new calibration, making use of available J—K 
colours, and a newly-trained artificial neural network procedure 
(see Snider et al. 2001), to provide an alternative initial esti- 
mate of [Fe/H]. This latter step was taken because the HERES 
targets include objects that are known to be carbon-enhanced; 
the strong CH G-bands in these objects could possibly perturb 
the B — V colours used in the Beers et al. calibration. We take 
the average of the two determinations as our initial guess for 
[Fe/H]. The differences between the adopted initial guess and 
the final [Fe/H] values from the automated analysis of the high- 
resolution spectra are plotted in Fig. |5J and application of ro- 
bust methods (see, e.g., Beers, Flynn, & Gebhardt 1990) to the 
set of differences in the initial and final metallicity determina- 
tions yield estimates of the mean offset and standard deviation 
of +0.04 dex and 0.27 dex, respectively. The mean estimated 
relative error for the high-resolution abundance determinations 
of [Fe/H] is 0.12 dex (see Sect. l3.4l and Table|2j), indicating that 
the medium-resolution estimates have relative errors of about 
0.24 dex. Note that the [Fe/H] estimates for the worst outliers 
were often based on rather noisy medium-resolution spectra. 

An initial estimate of the surface gravity, log g, is also re- 
quired, which is refined in the automated analysis. For stars 



Fig. 5. Difference between the final [Fe/H] as deter- 
mined from the spectrum and the initial guess deter- 
mined from medium-resolution follow-up spectroscopy, where 
A[Fe/H]=[Fe/H](initial)-[Fe/H](final), plotted against T cff 
and [Fe/H](final). 

on the subgiant and giant-branch stages of evolution, it is well 
known that log g correlates very well with effective tempera- 
ture. Hence, in order to obtain a first-pass estimate of surface 
gravity we used the reported temperatures and surface gravities 
by Honda et al. (20()4i to derive the regression relation: 

log g = -9.301 + 2.273 T cff /1000 (1) 

Note that the Honda et al. program considers stars over a very 
similar range of metallicities and temperatures as the HERES 
program. For warmer stars, this fit leads to overly high initial 
estimates, so we limited the initial guess to a maximum value of 
log g = 4. That is, if T c g > 5850 K an initial guess of log g = 4 
is adopted, which was found to be adequate. Figure [6] shows 
the differences between the adopted initial estimate and the de- 
termination obtained from automated analysis of the snapshot 
spectra. 

As will be described in the next section, a ID analysis is 
employed and turbulent line broadening is modelled via the 
classical microturbulence £ and macroturbulence u ma cro pa- 
rameters. These parameters are derived from the spectrum. 
Initial guesses of £ = 1.8 km s _1 and u ma cro = 1-5 km s _1 
were employed. The final values of £ derived from the spec- 
trum are given in Table [2] and plotted in Fig. showing the 
usual correlation with gravity. Note that the stars lying above 
the main trend, around log g ~ 2 are the red horizontal branch 
stars (see Fig.[0- 



6 



P. S. Barklem et al.: The Hamburg/ESO R-process Enhanced Star survey (HERES) 



2 r 1 

CO 1 r \ • ^ 

01 1 L •« • • • ^ 

< = = 

-2 I ' j 

4000 4500 5000 5500 6000 6500 7000 
T e ,, [K] 

Fig. 6. Difference between the final log g as determined from 
the spectrum and the initial guess as described in the text, 
where Alogg — logg(initial) — log 3 (final). 
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Fig. 7. Plot of final microturbulence values £ as derived from 
the spectrum against gravity. Note that the stars lying higher 
than the general trend are the red horizontal branch stars (see 
Fig.EJ. 

3. Automated Spectrum Analysis 

3. 1. Analysis code 

Software for automated analysis of the spectrum has been de- 
veloped, based on the Spectroscopy Made Easy (SME) package 
by Valenti & Piskunov ( 1996 1. SME consists of three compo- 
nents: a spectrum synthesis component written in C++, a pa- 
rameter optimisation component written in IDL, and a user in- 
terface written in IDL. In this work we make use of the first two 
components. Our developed software is written in IDL, and es- 
sentially provides an alternative interface to the parameter op- 
timisation component, which in turn calls the spectrum syn- 
thesis component. We also made some minor adaptations and 
improvements to the SME codes, the most important of which 
are discussed below. The new interface software takes our in- 
put data and provides it to the parameter optimisation com- 
ponent without any need for user interaction. The input data 
in this case consist of a pipeline-reduced count spectrum with 
corresponding measurement error spectrum (corrected to the 
stellar rest frame), initial guess stellar atmosphere parameters 
and abundances (see Sect. 13.31 . and a list of spectral features 



with corresponding atomic or molecular data (see Sect. 13 .21 . 
The main tasks of the new interface software are, in a com- 
pletely automated fashion, to extract relevant spectral regions 
appropriate for the line list, to identify continuum points and 
normalise the spectra relative to the continuum, to make small 
adjustments (within the error of the wavelength calibration) to 
line central wavelengths such that they match the observed line 
centres, and to reject lines polluted by artifacts such as cosmic 
ray hits. The continuum points are identified by an automated 
procedure where for each spectral region required for the anal- 
ysis, the continuum is determined for a wide region, including 
the desired spectral region and 7 A on each side by iteratively 
fitting a low order polynomial to the count spectrum and dis- 
carding points more than one standard deviation below the fit- 
ted line in the subsequent iteration, thereby converging to the 
continuum high points. Cosmic ray hits are identified by being 
significantly above this fit. 

Given the input data in the correct form, the SME parame- 
ter optimisation code can then solve for any desired model pa- 
rameters. The parameter optimisation code uses the Marquardt 
algorithm (Marquardt 1 1 9631 Press et al. 1992 1 to obtain esti- 
mates of the parameters, through minimising the \ 2 statistic 
comparing model and observed spectra. 

The spectrum synthesis assumes LTE and a ID plane- 
parallel model of the atmosphere, where turbulence is mod- 
elled through the classical microturbulence and macroturbu- 
lence parameters. Note that in recent versions of SME, molecu- 
lar line formation is supported (see Valenti et al. 1998 1, and for 
speed the radiative transfer is solved using a Feautrier scheme 
(e.g. Mihalas 119781 not a Runge-Kutta scheme (cf. Valenti 
& Piskunov 1996 1. The original spectrum synthesis code did 
not account for scattering in the source function, i.e. contin- 
uous scattering is treated as absorption and thus S v = B v . 
Rayleigh scattering by neutral hydrogen can be an important 
contributor to continuous opacity particularly in the ultravio- 
let for metal-poor stars, especially giants, and thus it is im- 
portant to treat continuous scattering correctly in the source 
function (e.g. Griffin et al. 1982 1 . It should be computed as 
S v = (1 — e u ) J v + e v B v where e„ is the photon thermalisa- 
tion probability given by n v J (k v + <j v ) where k v is the absorp- 
tion coefficient (including line absorption) and o v is the contin- 
uum scattering coefficient. This has been implemented into the 
SME spectrum synthesis code, where the mean intensity J u has 
been solved using the perturbation technique of Cannon J197U . 
Further minor improvements included correction of Thorium 
partition functions (see Paper I), and the ability to handle data 
for line broadening by neutral hydrogen collisions from the 
Anstee, Barklem & O'Mara theory (see Sect. 13. 2> . 

A grid of ID, plane -parallel, LTE models covering the stel- 
lar atmosphere parameter space of the stars of interest, namely 
low metallicity F, G, and K stars from the main sequence 
to the giant stage, was computed using the 1997 MARCS 
code (Asplund et al. 119971 Gustafsson et al. I1975> . All mod- 
els use scaled solar abundances with the exception of the 
alpha-elements, which are enhanced by 0.4 dex. Throughout 
this paper solar abundances means those of Grevesse & 
Sauval (1998 1, except C which is taken from Allende Prieto 
et al. (20021. Microturbulence in the line blanketing of 2 km 
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s _1 was adopted as most of our stars are giants (see Fig.0. 
This should not have a great effect in metal-poor stars in any 
case. This model grid was incorporated into SME, allowing 
SME to solve for atmospheric parameters, where specific mod- 
els are obtained by interpolation. 

In Sect. l3.3l we will describe how we apply this code. First 
we describe the adopted line list. 

3.2. Line list 

A list of spectral lines and corresponding atomic and molec- 
ular data, suitable for our sample of stars, observational data, 
and analysis technique is required. For our automated line pro- 
file analysis technique, spectral windows where the model and 
observed spectrum are to be compared need also to be defined. 

The lines and spectral windows used are listed in 
Tables lA~TllA.5l in Appendix|X] along with comments on each 
element. Note that in order to avoid the cores of strong lines, 
any wavelength points in the chosen spectral windows where 
the observed flux is less than 0.5 of the continuum are auto- 
matically rejected. The line list and windows were, for the most 
part, built by selecting suitable lines from previous studies, par- 
ticularly Mc Willia m et al . (I1995alll995bl . Norris et al. 0996 1, 
and Sneden et al. ( 1996). As our technique is automated, our 
aim is to build a line list which requires either no adjustment 
from star to star, or possibly a small amount of adjustment that 
can be automated based on criteria that are known in advance 
of the analysis, such as stellar temperature. We must therefore 
be more selective than in interactive analyses. One must keep 
in mind, however, the need to have as many lines as possi- 
ble to provide better statistics, which is particularly important 
for low-quality spectra. Ideally we would like to choose lines 
which/or all stars of the type we are examining are unblended, 
or which are partially blended in such a way that they can be 
used if an appropriate comparison window is chosen that ne- 
glects the blended part of the line (usually one wing of the line). 

Our first step for selecting the lines was to apply all can- 
didate lines to an automated abundance analysis of four stars 
viewed to be extreme cases of possible blending for metal- 
poor stars. We used HD 20, moderately r-process enhanced, 
which is one of the most metal-rich of our stars, non-carbon- 
enhanced; CS 31082-001, an r-process enhanced, non-carbon- 
enhanced star; CS 22892-052 an r-process enhanced, carbon- 
enhanced star; and HE 0338-3945, an r- and s-process en- 
hanced, carbon-enhanced star. Based on examination of the 
four spectra we chose appropriate spectral windows for each 
line, rejecting lines where a suitable window for all four stars 
could not be easily chosen. We then performed an abundance 
analysis of each star and rejected lines which were not well fit 
by the derived abundance with respect to the majority of lines 
of that element, indicating most likely an unresolved blend, or 
possibly poor atomic data. Special attention was paid to lines 
that were blended in the carbon-enhanced stars, but not in the 
carbon-normal stars. We chose to only make use of these lines 
for the carbon-normal stars, as described below. These lines 
are marked with an asterisk in the line list. To err on the side of 
caution, we also used some other very carbon-enhanced stars in 



the sample, which we deemed too carbon-rich for our present 
method, to indicate other lines which could be blended by car- 
bon features. Finally, the line list and windows were refined 
based on testing using the complete sample. 

We empirically determined that the lines noted above to be 
blended in the carbon-enhanced stars, could be safely applied 
in stars where the G band has a depth not greater than 0.6 of the 
continuum flux. Further, in the warmer stars some chosen lines 
lie in the wings of high Balmer lines. Thus, we empirically 
derived a scheme for application of the line list to a given star. 
The list is applied to all stars with the following adjustments: 

- if the maximum depth of the G band of CH has a depth 
greater than 0.6 of the continuum flux, the lines seen or 
suspected to be blended in carbon-enhanced stars (marked 
with asterisks in the tables) are removed from the list; 

- lines close to Balmer lines are removed depending on the 
stellar temperature. 

These adjustments are fully automated. 

While the employed spectral lines were chosen from the 
above mentioned previous studies, we compiled our own 
atomic data for these lines, though often guided by these 
works. In particular, we attempted to update the data. The 
adopted oscillator strengths are discussed for each element in 
Appendix El The VALD database 4 (Kupka et al. fl999b and 
NIST Atomic Spectra Database 5 (see also Wiese et al. 1969 
Wiese & Martin [19801 Martin et al. [jggg] Fuhr et al. [T988l 
were used extensively, particularly for wavelengths and excita- 
tion energies. Radiative broadening, Stark broadening and col- 
lisional broadening by neutral hydrogen are included, though 
the data have not been presented in the tables; they may be 
obtained from the authors on request. Radiative broadening is 
taken in all cases from VALD, which is supplied from the cal- 
culations by Kurucz ( 1995). Collisional broadening by neutral 
hydrogen is, where possible, described by the Anstee, Barklem 
and O'Mara theory (e.g. Anstee & O'Mara [79911 Barklem 
et al. 2000 and references therein). In the absence of such cal- 
culations the data are taken from VALD which again orig- 
inates from Kurucz (1995) using the van der Waals theory 
(Unsold [l955> with a detailed calculation of the long range in- 
teraction constant C§ . In the absence of data from both of these 
sources the classical van der Waals theory is used, where Cq 
is estimated by the usual Unsold approximation. In both the 
last two cases we apply an enhancement of a factor of 2. Stark 
broadening is taken from VALD where available, which often 
originates from Kurucz ( 1995 ). We note that since the lines are 
typically weak, collisional mechanisms are generally not im- 
portant, except in a few lines which attain significant strength 
(e.g. the strongest lines of Mg, Sr and Ba). 

Where possible we have considered the hyperfine structure 
of spectral lines. Hyperfine structure is typically seen in ele- 
ments with odd atomic numbers. Though Sr and Ba have even 
atomic numbers, in the solar system they have non-negligible 
abundances of isotopes with odd atomic masses, and thus hy- 
perfine structure must be considered in these elements also. 

4 http ://w ww. astro. uni vie . ac . at/~ vald/ 
http://physics.nist.gov/cgi-bin/AtData/main_asd 
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Isotopic mixtures and splitting have been considered in Sr, 
Ba and Eu. In all other elements the dominating natural iso- 
tope is assumed (or equivalently in the case of even atomic 
numbers that isotopes with even atomic masses dominate and 
they have negligible isotope shift). Details of adopted hyper- 
fine structure and isotopic ratios are discussed for each element 
in Appendix[A]where it is considered. Hyperfine structure and 
isotopic splitting have been omitted from the tables for com- 
pactness. They may be obtained on request from the authors. 

3.3. Procedure 

We now describe the procedure used for analysis of the spectra. 
The goal is to determine the required stellar atmosphere param- 
eters and abundances from the spectra in the best possible way. 
In Sect. 13.11 a code was described which can analyse a given 
spectrum to obtain best estimates of chosen model parameters. 
While the code can in principal analyse the entire set of spectral 
lines and obtain all desired stellar parameters and abundances 
simultaneously in one step, this is not done, since for cases with 
many free model parameters such a procedure is susceptible to 
errors if strong dual dependencies of the model spectra on pa- 
rameters exists. Thus, our procedure, like traditional analyses, 
separates the determination of stellar parameters using appro- 
priate spectral features from the determination of abundances 
as we now describe. 

First, the T c g from Sect. 12.11 is adopted and held fixed 
throughout. We experimented with refining T e g from the spec- 
tra (essentially excitation equilibrium of Fe I), but found that 
the spectra are not of sufficiently high quality given the spec- 
tral features available to us. Thus we need to solve for the 
atmosphere parameters surface gravity log g, microturbulence 
£, macroturbulence v macI o, and metallicity [Fe/H]. Once these 
parameters are known we then solve for the individual abun- 
dances. All stars are assumed to be slow rotators; we adopt 
v sin i = 1 km s _1 for all stars, and this has no effect on the fi- 
nal abundances. Since the resolving power of observations can 
vary from star to star, we fix the instrumental broadening at a 
value R — 50000, significantly higher than expected and allow 
the derived macroturbulence to compensate for the difference. 
Thus our derived v maa o includes both the macroturbulence and 
a degree of the instrumental profile. 

Our basic procedure is broken into three steps. In the first 
two steps we employ only the Fe and Ti lines of our list. In the 
first step initial estimates of T e g, log g, [Fe/H], £ and w maC ro as 
described in Sect. l2.1l and l2.2l are adopted. The initial guess Ti 
abundance is chosen to be a scaled solar value based on the 
[Fe/H] initial guess. We then solve for the Fe and Ti abun- 
dances, log £fo 6 and log eii to obtain better initial estimates 
of these parameters before solving for the remaining stellar pa- 
rameters. 

The second step repeats the first, but now logg, £ and 
u macro are also free parameters. The minimisation to solve 
for log g is in essence equivalent to performing an ionisation 
equilibrium procedure for Fe and Ti. The resultant parameters 

6 We define the abundance parameter for element X in the standard 
notation log ex = log(A r x/A r H) + 12 where Nx are number densities. 



(logs', fmacro, [Fe/H], log eTi) are then adopted and held 
fixed for stage three, the analysis of lines of elements other than 
Fe and Ti in the list, where the only free parameters are the ele- 
mental abundances. All the remaining elements are considered 
separately, as this is most computationally efficient. Solar abun- 
dances scaled by [Fe/H] are adopted as initial guesses, which 
was found usually to be satisfactory; on rare occasions where 
convergence was not achieved this could be solved by search- 
ing for an initial guess which resulted in convergence. Note that 
in all steps, the single abundance which gives the best global 
fit to the spectrum is derived. It is also worth noting that we 
adopt scaled solar abundances for all elements which we do 
not analyse, except O which we assume enhanced by +0.4 dex; 
these abundances enter the molecular equilibrium and continu- 
ous opacity calculations. 

An additional stage is performed for checking purposes. In 
this stage we analyse the Fe and Ti line list keeping all pa- 
rameters fixed to those derived in the above stages except the 
abundance which may now vary from line to line (in SME this 
can be achieved in practice by fixing the abundance and solv- 
ing for log gf for each line, and observing the change from 
the log <?/-value in the line list). This permits us to check for 
trends with excitation and ionisation stage, which might warn 
of deficiencies in the adopted model. 

Three further steps are then performed which are required 
for the error analysis, the details of which will be described 
below in Sect. 13 .41 and in Appendix 151 However, for the error 
analysis we will require partial derivatives describing the de- 
pendence of each abundance on Teg, log g and £. Thus we solve 
for all abundances with fixed stellar parameters except that T e g, 
log g and £ are offset by 100 K, 0.3 dex and 0.2 km s _1 in each 
step respectively. We note that these last three steps increase 
computing time significantly (they comprise more than 1/2 of 
the computing time). However, we view the error analysis as 
important enough to justify this increase. 

It is worthwhile to make some comment on computing 
times. With the large number of molecular lines, isotopic and 
hyperfine components, our line list results in approximately 
900 individual line components to be included in the spectrum 
synthesis, which is where the vast majority of computation 
is used. Model convergence typically requires approximately 
3.5A^+13 model spectrum calculations, where N is the number 
of free parameters (Valenti & Piskunov ll996l l. This means for 
the main abundance step (step 3) we require of the order of 15 
model spectrum calculations of the complete spectrum. This is 
repeated three additional times for the error calculations. With 
a modern workstation typical computation times per star are of 
the order of 2-3 hours. 

3.4. Elemental detections and errors 

A first important step before assigning an error to a derived 
abundance is to determine if there is in fact a reliable detection. 
A Gaussian is fit to each line and the approximate equivalent 
width W is determined. Following Norris et al. (2001 1, the er- 
ror in W is then computed as aw — \\fn / {R[S / N]) where 
R is the resolving power, n the number of pixels integrated to 
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obtain W, and S/N the signal to noise ratio per pixel. In all 
cases we adopted R = 20000 and adopted the full-width of 
the fitted Gaussian measured in pixels for n. By then comput- 
ing m = WJ <Jw we classify a ma detection of the line. For 
each element we count the number of lines detected at the So- 
level or better. If the element has at least one line detected at 
the 3<7 level we assign the abundance as a detection. Otherwise 
the element is regarded as undetected, and the derived abun- 
dance discarded. Due to the possibility of weak blending, even 
in high quality spectra, it was further required that a line must 
be deeper than 10% of the continuum flux for a valid detection. 
Due to decreased likelihood of blending in the redder parts of 
the spectrum, we relaxed this condition to 5% for lines red- 
der than 4500 A. A similar detection scheme was used for CH 
bands, though W and n for the considered band must be esti- 
mated by a different technique. 

As one of the main goals of this work is to look at abun- 
dance scatter at low metallicities, and given the relatively low 
quality of the data, it is vital to have error estimates which are 
as accurate as possible to distinguish real scatter from that due 
to uncertainties. Thus, we have put significant effort into ob- 
taining realistic error estimates. 

Usually in analyses of stellar spectra, systematic errors are 
expected to dominate. This is not necessarily the case here. 
For elements with a small number of spectral features, the low 
S/N means that random measurement errors can be signifi- 
cant. However, for elements with a significant number of spec- 
tral features, e.g. Fe, measurement errors are in fact quite small, 
and systematic errors will dominate. 

Our method for calculating errors in abundances and abun- 
dance ratios is detailed in Appendix 151 First, we should em- 
phasise that we make a distinction between absolute and rela- 
tive errors. In our discussion we will refer to the uncertainty in 
the absolute abundance or abundance ratio as the absolute er- 
ror, and the uncertainty in the relative abundance between stars, 
that is from star to star, as the relative error. Appendix [Bj pro- 
vides a formalism for calculation of both absolute and relative 
errors; however, throughout this paper, unless otherwise stated 
we discuss the relative error. This is of most interest as we are 
interested in distributions of abundances. The main difference 
between absolute errors and relative errors is that some sys- 
tematic sources of error are presumed to cancel in the latter, in 
particular, errors in oscillator strengths and some part of the er- 
rors incurred from modelling uncertainties such as LTE and ID. 
We note however, that as our sample covers such a large stellar 
parameter space, the cancellation of errors from the modelling 
will be only partial. 

The error estimates include contributions from uncertain- 
ties in the observations, atomic and molecular data, continuum 
fitting, stellar parameters, model atmospheres and inherent as- 
sumptions in the spectrum modelling such as LTE and ID. 
We have included as many sources of error as possible, how- 
ever, some of these uncertainties are difficult to estimate quan- 
titatively, particularly those related to the model atmospheres 
and modelling assumptions. Our abundances and error esti- 
mates are provided with the explicit understanding that they 
are based on traditional ID LTE models and spectrum synthe- 
sis. Corrections for 3D and non-LTE effects often should be 



applied by the user. For example, it is known that the reso- 
nance line of Al is subject to strong deviations from LTE (e.g. 
Baumiiller & Gehren 1997 1. While corrections for cases such 
as Al where only a single line is employed are relatively straight 
forward, we point out that for other elements one needs to con- 
sider the correction to the best fit to all spectral features used 
here. 

3.5. Tests of the automated method 

Before presenting the results for the sample, we present the re- 
sults of some tests on our automated method. In particular, we 
compare our results where possible with the literature, and per- 
form test calculations to assess the robustness of the automated 
method at low S/N. 

3.5.1. Comparison stars 

Five stars, namely HD 20, HD 221170, CS 22186-025, 
CS 22892-052 and CS 31082-001, which have been studied 
in detail with better quality observational material by others, 
were included in our pilot sample. Four of these were in- 
cluded in the sample for the specific purpose of use as com- 
parison stars, while CS 22186-025 has been recently observed 
as part of the ESO First Stars programme (Cayrel et al. 2004 1. 
We now present a comparison of our results with some of 
the literature studies. For CS 22186-025, CS 22892-052, and 
CS 31082-001, we compare with the studies of Cayrel et al., 
Sneden et al. J2003i . and Hill et al. J20021 . respectively. For 
HD 20 and HD 221 170 there are a number of different studies 
which have included these stars. For the purpose of clarity, for 
the abundances of these two stars we chose to concentrate our 
comparison on just one study, that of Burris et al. ( 2000 1, which 
includes both these stars and has reasonable overlap with our 
work in terms of the elements analysed. 

First, in Table[5]we compare our stellar parameters for these 
stars with those used in the literature. Our temperatures are seen 
to be systematically warmer than those found in the literature 
for the comparison stars by of order 100 K. As mentioned in 
Sect. 12.11 and has been noted by other authors for individual 
stars, (e.g. Hill et al. 120021 Sneden et al. l2003t . the redden- 
ing maps of Schlegel et al. (1998 1 give larger reddenings than 
those of Burstein & Heiles (1982). As shown in Fig. [3] this 
leads to systematic differences in T e g, of around 100 K. In 
Table [3] we also quote our temperatures using the lower red- 
denings of Burstein & Heiles. These temperatures agree well 
with the literature within error (of order ±100 K). The remain- 
ing parameters agree reasonably well, given the typical abso- 
lute error bars of <J\ OEg ~ 0.35 dex, crpe/H] ~ 0.10 dex, and 
(T£ ~ 0.15 km s _1 . Considering the effect of the systematic 
offset in T e ff we see, in line with expectations, that it usually 
leads to slightly higher log g values (+100 K — > +0.2 dex), and 
slightly higher metallicities (+100 K — > +0.1 dex). 

To demonstrate the quality of the stellar parameter solu- 
tions, Fig. [8] plots abundances derived for individual Fe and 
Ti lines (the checking step of the procedure) against lower 
level excitation potential and wavelength of the line. Lines of 
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Table 3. Comparison of our final stellar parameters with those adopted in the literature for the comparison stars. T e g derived 
with the reddening maps of Burstein & Heiles ( 1982 1 instead of Schlegel et al. ( 1998! are also given for comparison. 



Star 




This Work 






Literature 




Reference/note 






logs 


[Fe/H] 




T e ff 


logs 


[Fe/H] 








[K] 


[cm s~ 2 ] 




[kms" 1 ] 


[K] 


[cm s~ 2 ] 




[kms" 1 ] 




HD20 


5445 
5416 


2.39 


-1.58 


2.30 


5351 
5475 
5375 


2.8 
2.41 


-1.4 
-1.4 


2.0 
1.90 


Alonso et al. 1 1999 1 
Burris et al. 1 2000 1 
Fulbright & Johnson <2003l 
Ilsinp Burstein & Heiles I 1982 1 


HD 221170 


4648 
4548 


1.57 


-2.14 


2.22 


4410 
4425 
4500 
4475 


1.0 

0.9 
1.15 


-2.0 
-2.1 
-2.2 


1.5 

2.75 
2.40 


Alonso et al. fl999l 
Rurris et al 1 2000 1 
Fulbright <2000l 
Fulbright & Johnson <2003l 
Using Burstein & Heiles 1 1982 1 


CS 22186-025 


4985 
4887 


1.70 


-2.87 


2.14 


4900 


1.5 


-3.0 


2.0 


Cayrel et al. <2004l 

Using Burstein & Heiles 119821 


CS 22892-052 


4884 
4843 


1.81 


-2.95 


1.67 


4800 


1.5 


-3.1 


1.95 


Sneden et al. <2003l 

Using Burstein & Heiles 119821 


CS 31082-001 


4922 
4874 


1.90 


-2.78 


1.88 


4825 


1.5 


-2.9 


1.8 


Hill et al. (ME\ 

Using Burstein & Heiles 1 1982 1 



the neutral and singly ionised species are distinguished in the 
plots, and thus it is clearly seen that the derived log g values 
satisfy ionisation equilibrium for both Fe and Ti. The plots 
against excitation potential, particularly for Fe, demonstrate 
that our spectra are not of sufficiently high quality to specify 
T e f{ with the required precision as seen from the large scat- 
ter in the plot for CS 22186-025 which has S/N ~ 32. The 
plots show our temperatures are reasonably consistent with ex- 
citation equilibrium (noting of course that one may question 
the validity of ID LTE excitation equilibria); however, as dis- 
cussed in Sect. l2.ll we find that on average excitation equilibria 
would give a cooler temperature scale. Finally, the plots against 
line wavelength indicate no significant systematic effects in 
the bluer lines due to blending or continuum placement. We 
also checked for significant trends with line strength and found 
none. Similar results were obtained for all stars in the sample. 

Figure [9] compares abundances for our comparison stars 
with results from the literature. The results show generally rea- 
sonable agreement within error, thus supporting not only our 
method for obtaining abundances, but our error estimates as 
well. For CS 22892-052 and CS 31082-001 there is a system- 
atic offset which is a result of our slightly warmer tempera- 
tures. This is clearly demonstrated by Fig.[H)] where we com- 
pare results for CS 22892-052 and CS 31082-001 where we 
have reanalysed the spectra adopting the temperatures used in 
the relevant literature work. The remaining small differences 
are likely attributable to use of different lines and atomic data, 
particularly oscillator strengths and hyperfine structure. 

The differences for HD 20 and HD 221170 are not im- 
proved by adopting temperatures used by Burris et al. (|2000 1, 



1 o - HD 20 (Burris et al. 2000) 



: mean = -0.08, a = 0.12 



which originate from Pilachowski et al. (1996 1. For HD 20, 
our temperature is already in reasonable agreement with that of 
Burris et al., yet we obtain systematically lower abundances by 
of the order of 0.3 dex. We reanalysed the spectrum adopting 
the stellar parameters of Burris et al., and the new comparison 
shown in Fig. shows a significant improvement. Thus, the 
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Fig. 11. Comparison of abundances as in Fig. M for HD 20 
where we have adopted the same temperature, gravity and mi- 
croturbulence as Burris et al. (2003. 



main reasons for the differences are our significantly higher 
microturbulence and lower gravity. We found no evidence of 
a significant trend of Fe abundance with line strength or ion- 
isation stage to indicate our microturbulence or gravity solu- 
tions are in error. In spite of the fact that our stellar param- 
eters for HD 221170 differ significantly from those adopted 
in Burris et al., the abundances are in reasonable agreement, 
probably due to compensation of the T c s , log g and microtur- 
bulence differences. Note that this star has a large colour excess 
E(B — V) « 0.12 in our sample, and a large discrepancy with 
the Burstein & Heiles maps (see Fig. [3}. 

While these systematic offsets, of order 0.1-0.2 dex, are 
present and attributed to differences in stellar parameters, per- 
haps the most important result is the scatter in the difference, 
as this indicates the accuracy with which the abundance pat- 
tern in a given star is reproduced. We find for all comparisons 
scatter of 0.1 to 0.15 dex. Thus we conclude that the relative 
abundance patterns are accurate to of order 0.15 dex, within 
our quoted error bars. 
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Fig. 8. Plots showing trends of abundances found from individual Fe and Ti lines with excitation potential x (left column) and 
wavelength A (right column) for two comparison stars, CS 22186-025 S/N ~ 32 and CS 31082-001 S/N ~ 90. Circles and 
crosses indicate lines from the neutral and singly ionised species respectively. The dashed line in each case shows the determined 
abundance from the simultaneous best fit to all lines. 



3.5.2. Robustness at low S/N 

Since most of our comparison stars are relatively bright, the 
spectra are usually of higher S/N than is typical for this 
project. The lowest quality comparison star spectrum is that of 
CS 22186-025, which has on average S/N w 32 per pixel. As 
some of our spectra are of even lower quality, it is important 
to ascertain what effect S/N has on our results. It should be 
pointed out that, provided our detection classification method 
and error estimates are reasonable, the pertinent issues are: (a) 
what is the minimum precision in abundances that is useful, 
and (b) what S/N is necessary to achieve it. 

To investigate this we took the spectrum of CS 31082-001, 
which has S/N w 90, and degraded it to lower S/N values 
by multiplying with appropriate Poisson noise and rescaling so 
that count numbers are consistent with the noise level. These 
spectra were then run through the entire automated spectrum 
analysis process completely independently from the original 
spectrum. The differences in derived atmosphere parameters 
(noting the model [Fe/H] is determined by the Fe abundance) 
and abundances for four degraded spectra compared to the orig- 
inal spectrum are given in Tabled along with the relative error 
bar for that derived parameter. First, we note that almost every 
difference reported was less than our computed error bar for the 
degraded spectrum, giving confidence in our estimated errors. 
Secondly, the differences are generally reasonable for most el- 
ements, less than 0. 1 dex. However, it is immediately obvious 
that for most elements there is a systematic trend towards slight 
underestimation of the abundances. We attribute this system- 



atic trend to systematically lower continuum placement, which 
is to be expected, since at low S/N it is impossible to distin- 
guish weak lines from continuum. Therefore, there should be a 
natural trend towards the continuum being placed too low, and 
thus underestimation of abundances, as the weak lines make 
the continuum appear slightly lower. This problem could per- 
haps be circumvented if one a priori defined continuum regions 
from higher quality spectra to be used for normalising the spec- 
tra; however such a procedure would be dependent on the simi- 
larity of the spectra. As noted, the error in the abundance at this 
S/N is much larger in any case. 

Based on these results, we conclude that our method seems 
to be reasonably robust even at S/N as low as 15 for elements 
with a large number of lines, such as Fe and Ti. However, as one 
expects, detections for a number of elements cannot be made 
at such low S/N for stars of this metallicity. Elements with a 
small number of lines which are still detectable at low S/N at 
this metallicity, for example Mg, Ca and Ni, are less robust, 
but at S/N > 15 the induced error usually does not appear to 
exceed of order 0.1 dex. 

3.5.3. Example fits to spectra 

Examples of portions of the spectra, and fits to the data ob- 
tained, are shown in Fig. ^1 f° r an r-process enhanced star 
HE 1127-1143. The example was chosen to be representative 
of the sample; it is a giant star with S/N ~ 49, typical for the 
sample. 
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Fig. 9. Comparison of abundances for stars previously studied in the literature. Differences are in the sense of Aloge = 
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Fig. 10. Comparison of abundances as in Fig. [9] where we have adopted the same temperature as the literature study. 

4. Results 



The derived stellar parameters and elemental abundances for 
the sample are presented in Table |2] 7 Below, we note and dis- 
cuss objects of particular interest found in the sample. In the 
next section we will discuss the more general behaviour of the 
abundances and implications for understanding the chemical 
evolution of the Galactic halo. 

It is worth emphasising that some care must be taken in in- 
terpreting the results presented here involving elements which 
are not always, or even usually, detected in our spectra, which 
may lead to selection effects. Fourteen elements, C, Mg, Al, 



7 Table|2|in its entirety is available only electronically 



Ca, Sc, Ti, Cr, Mn, Fe, Co, Ni, Sr, Y and Ba, are almost always 
detectable at the 3a level in the spectra. Another eight elements 
are analysed, V, Zn, Zr, La, Ce, Nd, Sm and Eu, which can usu- 
ally only be detected in the spectra of the least metal-poor stars 
of our sample, spectra with higher than usual S/N, if the abun- 
dance is enhanced, or a combination of these factors. Thus, par- 
ticular care must be taken in interpreting results involving these 
eight elements as the data are incomplete. 



4. 1. Objects of Interest 

The complete sample consists of 253 stars, of which four 
were comparison stars already known to be r-process enhanced 
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Table 4. Comparison of results for CS 31082-001 using spectra with degraded S/N. The difference from our results using the 
original spectrum S/N 90 in Table [2] is quoted, with the relative error in the derived abundance for the degraded S/N in 
parentheses. A blank means a non-detection (below 3a) at this S/N for this star. 
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A Parameter 
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log£(Mg) 
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log£(Ca) 
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metal-poor stars. Among the remainder we have identified a 
number of interesting objects, which we will now summarise. 
For the discussion in this subsection (§ 14. li we consider only 
the 249 other stars as "the sample". Though CS 22186-025 has 
been recently observed by Cayrel et al. (2004) it is included in 
our statistics as it had not been observed at the time of selec- 
tion. 

First, this work adds significantly to the number of con- 
firmed very metal-poor stars. The [Fe/H] distribution for the 
sample is shown in Fig. [2] The sample contains 49 stars with 
[Fe/H] < -3, and 181 stars with -3 < [Fe/H] < -2. Only 
12 stars are currently known with [Fe/H] < —3.5 (Beers & 
Christlieb 2005 1 and we find just one new star in this regime, 
HE 1300+0157 with [Fe/H] = -3.76 ± 0.18. This star is at the 
base of the giant branch (T c s ~ 5400 K, log g ~ 3.4) and is 
carbon enhanced with [C/Fe] ~ 1.2. 

The main goal of the HERES survey is to identify r-process 
enhanced stars, to address questions about the r-process. We 
were able to detect both Eu and Ba at the 3cr level in 62 stars in 
the sample. Of these, 57 are judged to be "pure" r-process stars, 
as they have [Ba/Eu] < 0, see Fig.[0] In Paper I we made the 
distinction between r-I and r-II stars, 0.3 < [Eu/Fe] < 1.0 and 
[Eu/Fe] > 1.0 respectively (and [Ba/Eu] < in both cases), 
on the basis that previous work had suggested a bimodal dis- 
tribution of [Eu/Fe] with a lack of stars found in the range be- 
tween 1.0 and 1.5. A histogram of detected [Eu/Fe], is plotted 
in Fig. [2J however, due to the significant incompleteness of 
the Eu abundances and the detection bias towards high [Eu/Fe], 
only the high [Eu/Fe] side of this distribution is reliable, and 
thus this plot gives little information about the cosmic distri- 
bution. We will return to the question of the distribution of r- 
process enhancement in Sect. 15. 31 Following this classification 
system, of these 57 stars, eight are r-II stars ([Eu/Fe] > 1.0) 
which are listed in Table [5] (two of which were announced in 



Paper I), and 35 are r-I stars (0.3 < [Eu/Fe] < 1.0), while 
14 do not have r-process enhancement ([Eu/Fe] < 0.3). This 
implies a frequency of around 3% and > 14% for r-II and r- 
I stars respectively among metal-poor stars ([Fc/H] < —1.5). 
It must be borne in mind that these frequencies are dependent 
on detection in our spectra, and thus are certainly lower limits, 
particularly for r-I stars. For example, at the lowest metallici- 
ties, an r-I star would often not be detectable with our typical 
spectra. 

The abundance patterns for the neutron-capture elements 
for the eight r-II stars are shown in Fig. \^\ along with the 
scaled solar system r-process abundances, demonstrating that 
these elements follow this pattern in these stars. In fact, all 
57 pure r-process stars, as judged from Ba/Eu, follow the 
pattern quite well for Ba and heavier elements. The lighter 
neutron-capture elements, Sr, Y and Zr, tend to deviate from 
the solar-r-process pattern in stars of lower r-process enrich- 
ment, which will be discussed in Sect. 15.41 The scaled solar 
system s- and r-process abundances used throughout this pa- 
per are based on the s-process fractions from the stellar model 
of Arlandini et al. (1999 1, the assumption that the remainder 
of neutron-capture elements are produced by the r-process, 
and meteoritic abundances from Grevesse & Sauval (1998 1. 
We chose the Arlandini et al. (1999 1 fractions over those of 
Burris et al. ( 2000 1 for our comparisons, since we found consis- 
tently better agreement with Y abundances in the pure r-process 
stars. The abundance patterns are otherwise very similar for the 
other neutron-capture elements considered here. 

The eight r-II stars are all giants with T c s ~ 5100 K, and 
have a quite narrow range in metallicity, —3.2 < [Fe/H] < 
—2.6. The r-II stars, including the two previously known exam- 
ples CS 22892-052 and CS 31082-001, are centred on a metal- 
licity of [Fc/H] = —2.81, with a very small scatter, on the or- 
der of 0.16 dex. Note, we find one r-II star with [Fe/H] < — 3 
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Fig. 12. Examples of the spectra and fits for HE 1 127-1 143 (an r-process rich star), which has average S/N ~ 49, T e s = 5224 K, 
log <7 ~ 2.6, [Fc/H] ~ —2.7, and is thus reasonably typical for the sample. The double line is the observed spectrum and the 
single line the fitted model spectrum (noting that it is not only the plotted observations which are being fit, but all considered 
lines of the element). The shaded regions are the comparison windows used for fitting. 



Table 5. New r-II and s-II objects. Quoted error estimates are relative errors. N3 is defined in the caption to table[2] 
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r-II stars 
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5103 


2.45 


-2.81 ±0.13 


46 


0.18 ±0.17 


1 


0.34 ±0.19 




1.06 ±0.15 


2 


-0.71 ±0.17 


CS 29497-004 


5013 


2.23 


-2.81 ±0.13 


42 


0.22 ± 0.18 


1 


1.21 ± 0.15 




1.62 ± 0.15 


4 


-0.41 ± 0.17 


HE 0430-4901 


5296 


3.12 


-2.72 ±0.12 


46 


0.09 ±0.18 


1 


0.50 ±0.20 




1.16 ±0.17 


3 


-0.65 ±0.17 


HE 0432-0923 
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2.64 


-3.19 ±0.13 


41 


0.24 ±0.17 


1 


0.72 ±0.18 




1.25 ±0.15 


2 


-0.53 ±0.16 


HE 1127-1143 
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-2.73 ±0.14 


46 


0.54 ±0.17 
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0.63 ±0.19 




1.08 ±0.15 
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0.35 ±0.17 


1 


0.59 ±0.18 




1.05 ±0.15 
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s-II stars 
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-2.71 ±0.11 


22 


2.45 ±0.18 
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2 


1.81 ± 0.16 
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among 49 stars in our sample, a frequency of 2%, and 7 r-II 
stars among 118 stars in the range —3 < [Fc/H] < —2.5, 
a frequency of 6%. This perhaps indicates that r-II stars are 
more rare at [Fc/H] < —3; further study in a larger sample 
with better spectra would be desirable to resolve this question 
more definitively. All eight r-II stars have quite normal C/Fe 



abundances, the largest being [C/Fe] ~ 0.5, noting of course 
that the present sample is biased towards CH weak stars. The 
low C abundances mean that these stars are all candidates for 
studies of the actinide elements, particularly Th and U, in such 
stars. CS 22892-052 remains a unique object as the only r-II 
star known with [C/Fe] ~ 1. 
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Fig. 15. Abundance patterns for neutron-capture elements in the r-II stars of the sample. The full line shows the scaled solar 
system r-process abundances based on the r-process fractions from Arlandini et al. ( 1999), meteoritic abundances from Grevesse 
& Sauval ( 1998 1, scaled to match our observed Eu abundance. Plotted error estimates are relative errors. 



The majority of the r-I stars are also giants, however, four 
such stars are unevolved (i.e. not yet on the giant branch T e g > 
5500 K, log 3 > 3), namely HE 0341-4024, HE 0534-4615, 
HE 0538-4515 and HE 2301-4024. As seen in Fig.[l7]the r-I 
stars span the range —3.4 < [Fe/H] < —1.5, practically the 
entire metallicity range of our sample. We note that the four 
unevolved stars all have [Fc/H] > —2.1, but that there is a 
selection effect at work here; the Eu II lines observed will be 
weaker in unevolved stars than in evolved stars with the same 
abundances, and thus our survey preferentially detects Eu in 
giant stars. Our survey is also biased towards cool giants as 
described in Paper I. 

Five stars in which both Eu and Ba were detected are found 
to be s-process rich as judged from their having [Ba/Eu] > 
0.5. 8 All five are carbon-enhanced, and three are strongly Eu 
enhanced with [Eu/Fe] > 1.0. The three stars with [Ba/Eu] > 



8 As described in appendix|A| we have assumed r-process isotopic 
composition of Ba in all stars which may lead to underestimation of 
the Ba abundance and thus Ba/Eu in stars with significant s-process 
contributions. Similarly, we assume no 13 C which may lead to errors 
in the C abundances for such stars. 



0.5 and [Eu/Fe] > 1.0, are listed in Table|5] A number of sim- 
ilar stars are known; e.g. Hill et al. (2000), Aoki et al. ( 2002b I 
and Cohen et al. 



(20031. We shall refer to these stars as "s- 
II" stars, though we note that whether the neutron-capture el- 
ements in such stars are produced predominantly by the s- 
process or by both the s- and r-processes is a matter of cur- 
rent debate (e.g. Johnson & Bolte 2004 1. We note that these 
three stars all have large C enhancements [C/Fe] > 1.5, and 
are unevolved. It must be borne in mind, however, that since 
this sample is limited to the stars with spectra only weakly pol- 
luted by CH features, we were only able to analyse the warmest 
strongly C enhanced stars. The complete HERES sample con- 
tains 72 CH strong stars which were not analysed here, among 
which there will be a large number of s-process rich stars. As 
we will discuss in the next section, a further such s-II star is 
suspected based on other abundances, but is not confirmed as 
Eu is detected at below the 3a level. 

The remaining two s-process rich stars with [Ba/Eu] > 0.0 
yet more normal Eu enhancement ([Eu/Fe] < 0.6) have abun- 
dance patterns of neutron-capture elements which are reason- 
ably well matched by a scaled solar s-process abundance pat- 
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Fig. 14. Histogram of [Eu/Fe] values. The dashed line includes 
all stars where Eu was detected, including the four comparison 
stars. The full line includes the pure r-process stars only (i.e. 
[Ba/Eu] < noting Ba must be detected; for a small number 
of stars where Eu was detected, Ba was not detected due to 
cosmic ray pollution of the Ba II resonance line) and excludes 
the comparison stars. 



Fig. 13. A plot of [Ba/Eu] vs [Fe/H] (upper panel) and [C/H] 
(lower panel). The horizontal full line shows the pure solar r- 
process value, computed from the solar r-process fractions of 
Arlandini et al. ( 1999), and the dashed line the solar s-process 
value. The pure r-II stars are shown as stars, the r-I stars as 
diamonds, the pure r-process stars without excess r-process el- 
ements as triangles, the s-II stars as asterisks, and the two re- 
maining s-process-rich stars as squares. The average relative 
error bar is shown in the bottom left. 

tern. Figure^] shows the abundance patterns for the three s-II 
stars. In these cases the Ba, La, Ce and Nd abundances can 
be well reproduced by a scaled solar-system s-process abun- 
dance pattern. The Sr and Y abundances, however, are much 
lower. We notice, however, that these abundances are consis- 
tent with the scaled solar system r-process abundance pattern 
normalised to Eu. On the other hand, the s-process is expected 
to be skewed towards heavier elements in metal-poor environ- 
ments (e.g. Busso et al. 1999). A summary of suggested pos- 
sible production scenarios for the s-process rich metal-poor 
stars has been recently given by Johnson & Bolte (2004i. As 
stated by those authors, to distinguish different scenarios more 
strongly, measurements of as many heavy elements as possible 
are needed, and studies of these stars using better observational 
material will be the subject of future work as part of the HERES 
survey. 

5. Discussion 

5. 1. Abundance Trends and Scatter with Metallicity 

We now discuss the general behaviour of the abundances, par- 
ticularly trends and scatter in abundance ratios of interest. The 



abundance ratio trends of [X/Fe] with [Fe/H] and [X/Mg] with 
[Mg/H] have been plotted in Figs.ll7landll8lrespectivelv for el- 
ements with a reasonable number of detections. The estimated 
la scatter in the y-variable is shown, following the definition 
for the scatter used by Karlsson & Gustafsson ( 2005 ), such that 
at a given x-coordinate 32% of the stars lie outside the lcr scat- 
ter lines, 16% above and below. In practice, this has been com- 
puted for a given x-coordinate by taking the n nearest stars in 
the x-variable, where n is the lesser of n — 50 or n stars /3, 
with n stars being the number of stars in the diagram. The re- 
sults are smoothed over the average width of bins containing 
n stars to remove transient behaviour caused by outliers. In 
Table [6] we compare the mean measured scatter in [X/Fe] and 
[X/Mg], across the range of [Fe/H] and [Mg/H] respectively, 
with the average error estimates. The minimum and maximum 
measured scatters are also reported to give an indication of the 
range of variation with [Fe/H] and [Mg/H] . 

For Mg, Ca, Sc, Ti, Cr, Fe, Co, and Ni, the measured scat- 
ters are slightly less than the scatters expected from the typ- 
ical estimated relative error in the abundance ratios, suggest- 
ing that our computed errors overestimate the real relative er- 
ror. Bearing in mind the difficulties in accurately calculating 
the errors, this suggests that the cosmic scatter is small (or 
even non-existent) at these metallicities. These results are in 
line with previous results for smaller samples with smaller 
uncertainties for some of these elements, such as the stud- 
ies by Arnone et al. (|2005 1, Cayrel et al. J2004> and Cohen 
et al. (2004). It should be noted that in all these elements, ap- 
parent outliers may perhaps be explained simply by the errors 
in the analysis; in a sample of 253 stars outliers due to ran- 
dom errors are expected. In particular, examining the outliers 
we noted that they often showed trends in the abundances de- 
rived from individual Fe I lines with excitation that were at the 
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Fig. 16. Abundance patterns for neutron-capture elements in the s-process rich stars of the sample with strong Eu enhancement, 
[Eu/Fe] > 1.5, the s-II stars. The full line shows the scaled solar system r-process abundances based on the r-process fractions 
from Arlandini et al. ( 1999 1, meteoritic abundances from Grevesse & Sauval ( 1998 ), scaled to match our observed Eu abundance. 
The dashed line shows the equivalent scaled solar system s-process abundances, scaled to match our observed Ba abundance. 
Plotted error estimates are relative errors. 



edge of the distribution of such slopes for the entire sample 
(see Fig.|4j, and occur in elements where particularly tempera- 
ture sensitive spectral features have been employed, such as Co 
and Ni. Thus, we believe that more often than not for these ele- 
ments, the outliers are more likely due to errors in the analysis, 
such as T c ff being in error, than indicative of any real over- or 
under-abundance of the element in question. 

The results for Al, V, Mn and Zn are less clear. In the case 
of Al, the scatters are marginally larger than the typical error 
bars; however, we note the Al abundances are based on a sin- 
gle line and are very sensitive to the S/N of the spectrum (see 
Table[4}. Further, this line is known to be subject to deviations 
from LTE. Both these facts may lead to increased scatter in 
the derived abundances. In the cases of V and Zn, the number 
of detections is small and the abundances for these elements 
are based on quite weak features and are therefore suscepti- 
ble to overestimation due to unresolved blends. For Mn, the 
scatter at low metallicity does marginally exceed the error es- 
timates, with the hint of a "bump" of Mn enhanced stars at 
around [Fc/H] ~ —3. A similar feature is suggested by the 
results of Cayrel et al. J2004t . 

The elements C, Sr, Y, Ba and Eu, and perhaps Zr, show 
scatter at low metallicities, [Fe/H] <, -2.5 and [Mg/H] ^ 
—2.2, significantly larger than can be explained from the er- 
rors in the analysis, implying the scatter is cosmic in ori- 
gin. At higher metallicities the scatter among the non-C- 
rich stars, does not greatly exceed that expected from the 
errors in the analysis, but perhaps indicates some cosmic 
scatter. Note, the scatter in [Ba/Fe] and [Ba/Mg], and to 
a lesser degree Sr and Y, at higher [Fe/H] and [Mg/H] is 
affected by C-enhanced stars which have high Ba/Fe and 
Ba/Mg ratios, see Fig. [28] The scatter at low metallicity 
seems, in most cases, to be increasing monotonically with 
decreasing metallicity. Similar results have been found, for 
example, by McWilliam et al. ( U995b> . McWilliam (ll^98l . 
Norris et al. (I200H . and Bums et al. J2000> . The results of 
Norris et al. {2001 1 suggested the existence of a larger scatter 
in [Sr/Fe] than for [Ba/Fe] at the lowest [Fe/H]. We do not find 
this result; the scatters are of comparable magnitude. Note that 
in contrast to the results for lighter elements discussed above, 



the cases of under-abundances and over-abundances of these 
elements with respect to the general trends are sometimes sig- 
nificant. 

In summary, among the stars without strong C enhance- 
ment, at about [Fe/H] ^ —2.5 or equivalently [Mg/H] 
—2.2, our analysis indicates that the cosmic scatter in all abun- 
dance ratios is small. This implies that at around this level of 
enrichment the Galactic halo was reasonably well mixed. At 
lower metallicities C, Sr, Y, Ba and Eu all show evidence for 
real cosmic scatter, while the results for the remaining elements 
still admit only a small amount of cosmic scatter within the 
errors of our analysis. This dichotomy implies that while the 
lack of scatter in these elements at [Fe/H] ^ —2.5 might be 
explained by a well mixed Galactic halo, the small scatter at 
[Fe/H] ^ —2.5 has a different explanation. One possibilty is 
suggested by the stochastic models of metal-poor enrichment 
by Karlsson & Gustafsson (2005 1. The small scatter among the 
most metal-poor stars may be tentatively explained by cosmic 
selection effects in contributing supernova masses, and a rel- 
atively narrow range of masses of gas the newly synthesised 
elements mix with (so called mixing masses). Alternatively, as 
suggested by Arnone et al. (2005), the scatter in C, Sr, Y, Ba 
and Eu at these metallicities might be explained by the addition 
of sources of these elements which produce negligible amounts 
of the other elements such as Mg and Fe, such as low mass su- 
pernovae of type II. 

5.2. Carbon 

The work of Bromm & Loeb ( 2003 1 suggests that halo stars 
with low C abundances, [C/H] ~ —3.5 according to their 
model, and low oxygen abundances, are probably true second 
generation stars. In Fig. ^] we plot the [C/H] values for our 
sample, identifying evolved (on the giant branch or red hori- 
zontal branch) and unevolved stars (all others). Though we do 
not have oxygen abundances, we see the majority of stars lie 
above the —3.5 level. It would be interesting to obtain O abun- 
dances for the stars lying near or below [C /H] ~ —3.5 to iden- 
tify likely second generation stars among this sample. 
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Fig. 17. Abundances ratios [X/Fe] plotted against [Fe/H] for elements with significant numbers of detections. Full lines show 
estimated la scatter. The average relative error bars are shown in the bottom left. Note differing scales on the y-axes. 
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[Mg/H] [Mg/H] [Mg/H] 




Fig. 18. Abundances ratios [X/Mg] plotted against [Mg/H] for elements with significant numbers of detections. Full lines show 
estimated la scatter. The average relative error bars are shown in the bottom left. Note differing scales on the y-axes. 



20 



P. S. Barklem et al.: The Hamburg/ESO R-process Enhanced Star survey (HERES) 



Table 6. Comparison of measured scatter in Figs. 1171 and 1181 
with relative errors in the abundance ratios. For each plot we 
report the minimum, mean and maximum measured la scatter 
fmeas across the range of [Fe/H] or [Mg/H], which is com- 
pared to the average relative error in the abundance ratio er rcl . 
The ratio of the mean measured scatter to the estimated error 
(cmeas) /c rel is reported in the last column. 





min((T moas ) 


(^raeas) 


max(cr mcas ) 


(7 


(o"mcas)/0" rCl 


[C/Fe] 


0.22 


0.32 


0.53 


0.18 


1.73 


[Mg/Fe] 


0.07 


0.09 


0.12 


0.15 


0.59 


[Al/Fe] 


0.17 


0.20 


0.21 


0.17 


1.13 


[Ca/Fe] 


0.08 


0.09 


0.10 


0.15 


0.59 


FSr/Fel 


0.09 


0.11 


0.13 


0.16 


0.69 


[Ti/Fe] 


0.08 


0.09 


0.12 


0.15 


0.61 


rV/Fe] 


0.06 


0.09 


0.11 


0.15 


0.60 


[Cr/Fe] 


0.10 


0.12 


0.14 


0.16 


0.72 


[Mn/Fe] 


0.13 


0.16 


0.21 


0.16 


1.03 


[Co/Fe] 


0.10 


0.11 


0.13 


0.16 


0.71 


[Ni/Fe] 


0.12 


0.14 


0.17 


0.18 


0.80 


[Zn/Fe] 


0.08 


0.15 


0.19 


0.15 


1.01 


[Sr/Fe] 


0.19 


0.32 


0.59 


0.19 


1.71 


[Y/Fe] 


0.21 


0.25 


0.32 


0.17 


1.47 


[Zr/Fe] 


0.13 


0.19 


0.23 


0.16 


1.18 


[Ba/Fe] 


0.37 


0.50 


0.69 


0.18 


2.77 


[Eu/Fe] 


0.21 


0.37 


0.62 


0.16 


2.34 


[C/Mg] 


0.24 


0.35 


0.51 


0.18 


1.89 


[Al/Mg] 


0.15 


0.19 


0.23 


0.18 


1.03 


[Ca/Mg] 


0.07 


0.08 


0.09 


0.15 


0.56 


[Sc/Mg] 


0.10 


0.11 


0.12 


0.18 


0.60 


[Ti/Mg] 


0.08 


0.09 


0.11 


0.17 


0.55 


[V/Mg] 


0.12 


0.13 


0.14 


0.16 


0.79 


[Cr/Mg] 


0.12 


0.15 


0.19 


0.18 


0.84 


[Mn/Mg] 


0.17 


0.21 


0.30 


0.17 


1.26 


[Fe/Mg] 


0.08 


0.10 


0.14 


0.15 


0.63 


[Co/Mg] 


0.11 


0.13 


0.16 


0.17 


0.79 


[Ni/Mg] 


0.15 


0.17 


0.19 


0.19 


0.88 


[Zn/Mg] 


0.09 


0.18 


0.26 


0.15 


1.16 


[Sr/Mg] 


0.18 


0.29 


0.55 


0.21 


1.41 


[Y/Mg] 


0.23 


0.25 


0.30 


0.19 


1.36 


[Zr/Mg] 


0.15 


0.22 


0.29 


0.17 


1.27 


[Ba/Mg] 


0.33 


0.48 


0.60 


0.20 


2.40 


[Eu/Mg] 


0.22 


0.39 


0.60 


0.18 


2.18 



5.3. Heavy Neutron-Capture Elements, Ba-Eu 

Figure^] which we used earlier as a measure of r- vs s-process 
enrichment, plots [Ba/Eu] against Fe and C abundances. The 
plots show a clear separation between two groups in the halo, 
a separation which correlates with C enrichment. This distinc- 
tion was first seen in Mc William (1998), though with fewer 
stars. The scatter among the pure r-process stars is consistent 
with the observational uncertainties and we thus conclude that 
the cosmic scatter in Ba/Eu among pure r-process halo stars 
is small. This implies that the Ba/Eu abundances produced 
by the r-process in the early Galaxy are universal. Taking the 
weighted mean (weights are based on absolute errors), we find 
for the pure r-process stars analysed in this work ([Ba/Eu]) = 
—0.58 ± 0.03, where the quoted error is the weighted standard 



x -2 
o 

-3 




5 F , 



[C/Fe 



Fig. 19. A plot of [C/H] versus [C/Fe] following Bromm & 
Loeb 



J2003). Giants and red horizontal branch stars are plot- 
ted as open circles, while the subgiants and dwarfs are plotted 
as filled circles. 

deviation. The solar system r-process abundance ratio, based 
on the data from Arlandini et al. ( 1999 1, is —0.65. Our result 
seems to be in disagreement with that found at low metallici- 
ties by Truran et al. (2002 1, who compiled data from a number 
of studies, and found a large scatter in Ba/Eu. They noted the 
difficulties in analysing Ba in metal-poor spectra due to depen- 
dence on hyperfine and isotopic structure and micro turbulence. 
This is particularly problematic when combining data from dif- 
ferent studies. Our data on the other hand are homogeneously 
analysed, but the Eu abundances are incomplete and biased to- 
wards stars with strong r-process enhancement. 

It has been suggested (e.g. Burris et al. 120001 Truran 



et al. 2002 Simmerer et al. 2004) that La may be a better alter- 
native as a tracer of the s-process enrichment. Unfortunately, 
as the lines available to us are typically weak, La is often 
undetected in our low S/N spectra according to our detec- 
tion criteria. When plotting [La/Eu] vs [Fe/H], Fig.|20| we see 
no trend or significant scatter. Among the pure r-process stars 
([La/Eu]) = —0.40 ± 0.05. The solar system r-process abun- 
dance ratio is —0.38. 

As mentioned in Sect. 14.11 the distribution of r-process 
enhancement in the early Galaxy is of great interest. This is 
usually traced by [Eu/Fe], but our data are incomplete due 
to the inability to reliably detect low Eu abundances in our 
snapshot spectra. One may, however, try to reconstruct the 
r-process enhancement by instead using the Ba abundances, 
which are practically complete due to the ease of detecting the 
Ba II resonance line at 4554 A. Following Raiteri et al. {1999 1 
and Burris et al. 



(2000), where possible we attempt to recon 
struct the r-process-only contribution to Ba. We define the ra 
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Fig. 20. A plot of [La/Eu] vs [Fe/H]. Symbols and lines have 
the same meanings as in Fig. 1131 
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tio [Ba/Fe] r _ procoss = log(iV Bair /iV Fe )*-log(iVBa,r/iV F e) , 
where A^Ba.r indicates the r-process-only Ba abundance. Note, 
the zero point of the scale is set at the solar r-process-only Ba 
abundance, not the total solar Ba abundance. For the pure r- 
process stars, [Ba/Fe] r _ proccss = [Eu/Fe], since the Ba/Eu 
r-process ratio is fixed at the solar r-process-only ratio and Eu 
is produced practically entirely by the r-process. Previous stud- 
ies indicate that the main s-process does not become important 
until [Fe/H] > -2.5, e.g. Burris et al. and Truran et al. d2002l . 
Thus we assume that at lower metallicities all Ba in non- 
C-rich stars is produced by the r-process. That is, for stars 
with [C/Fe] < 0.3 and [Fe/H] < -2.5, [Ba/Fe] r _ proce88 = 
[Ba/Fe] + 0.721. Those stars that are not pure r-process stars 
with higher C enrichment and higher metallicity are disre- 
garded. 

The results for [Ba/Fe] r „ process are plotted against [Fe/H] 
in Fig. Note, since we can only estimate [Ba/Fe] r _ proccss 
in stars with [Fe/H] > —2.5 if Eu is detected, there is an ob- 
servational bias here; we will miss stars with low r-process en- 
hancements at [Fe/H] > —2.5. The results of Burris et al., 
which cover this metallicity regime and are much more com- 
plete, indicate the scatter in the region —2.5 < [Fe/H] < —1.5 
is larger than seen in our results. They find stars in this regime 
with [Eu/Fe] ranging from as low as ~ —0.4 (significantly 
lower than in our data) to as high as ~ 0.8 (similar to our data). 
In any case, even accounting for the fact that this bias leads to 
an underestimation of the scatter at [Fe/H] > —2.5, the scatter 
at [Fe/H] < —2.5 in our data is significantly larger than seen at 
higher metallicity in the results of Burris et al. If we consider 
only the upper envelope of the abundance distribution, which 
should be well defined by our sample, there is an rapid transi- 
tion from a maximum value of [Eu/Fe] of ~ 1.0 to ~ 1.7 at 
[Fe/H] ~ —2.5. Thus, Fig.l21land the occurrence of r-II stars 
only at [Fe/H] < —2.5, suggests the r-II stars are extreme cases 
of a wide range of r-process enrichment in the early, chemically 
inhomogeneous Galaxy. Note, the cutoff of [Eu/Fe] = 1.0 be- 
tween r-I and r-II is suggested by the maximum enrichment 
levels at higher metallicity. 



Fig. 21. A plot of [Ba/Fe] s vs [Fe/H]. Symbols have the 

same meanings as in Fig. \T^\ For the closed symbols (pure r- 
process stars) [Ba/Fe] r _ proccss = [Eu/Fe]. Open circles show 
stars where [C/Fe] < 0.3 and [Fe/H] < -2.5, for which we 
assume [Ba/Fe] r _ process = [Ba/Fe] + 0.721. Note the scatter 
is underestimated at [Fe/H] > —2.5; see text. 

5.4. Light Neutron-Capture Elements, Sr, YandZr 

The production of light neutron-capture elements (partic- 
ularly Sr, Y, Zr) versus heavy neutron-capture elements 
(such as Ba, Eu), has become a topic of interest due 
to evidence of production of the former without signif- 
icant production of the latter, e.g. McWilliam (1998 1, 
Burris et al. (2000l, Truran et al. {2002 1, Travaglio et al. (I20041 
and Aoki et al. ( 12005V Figures l22l and 1231 plot [Sr/Ba] and 
[Sr/Eu] against Fe and C abundances. Significant scatter is 
seen in Sr/Ba at low metallicity, as found by McWilliam and 
Burris et al. In both cases, Sr/Ba and Sr/Eu, even among the 
pure r-process stars, a significant amount of scatter is seen at 
[Fe/H] < —2.5. For Sr/Ba, the scatter appears to increase 
quite uniformly with decreasing [Fe/H], the data showing quite 
clear upper and lower boundaries, apart from a small number 
of outliers which are usually C enhanced. Sr/Eu shows sim- 
ilar tendencies when the s-process rich stars are disregarded, 
though not as clearly due to the smaller number of stars. The 
r-II stars all have similar Sr/Ba and Sr/Eu ratios, which are al- 
ways among lowest of the non-C-enhanced stars. For Sr/Ba the 
weighted mean of the results for r-II stars, is ([Sr/Ba]) r _n = 
-0.44 ± 0.08, and for Sr/Eu ([Sr/Eu]) r _n = -0.98 ± 0.09. 
The solar system r-process values are —0.24 and —0.89 respec- 
tively. 

Truran et al. (j2002 1 examined the variation of Sr/Ba with 
r-process enrichment, as traced by Ba/Fe and Eu/Fe. Figure l24l 
shows Sr/Ba with Ba/Fe, and Fig. 125 1 shows Sr/Eu with Eu/Fe. 
We find similar results; however, with our large and homoge- 
neously analysed sample the scatter is well defined, and we see 
a trend for decreasing scatter in Sr/Ba and Sr/Eu with increas- 
ing r-process enrichment. 

Similar plots are shown for Y/Ba in Figs. 1261 and 1271 The 
results are similar to those for Sr/Ba, in particular, increas- 
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Fig. 22. A plot of [Sr/Ba] vs [Fe/H] (upper panel) and [C/H] 
(lower panel). Symbols and lines have the same meanings as in 
Fig- El wim circles representing stars where Eu is undetected. 

ing scatter in Y/Ba with decreasing metallicity and decreasing 
heavy r-process enrichment, and similar Y/Ba among the r-II 
stars. For the r-II stars we find ([Y/Ba]) r _ n = -0.47 ± 0.08; 
the solar system r-process value is [Y/Ba] = —0.48. 

Figures|22]|24]|26]and[^]identify two stars of interest, two 
metal-poor stars with [Sr, Y/Ba] —1 which are not carbon- 
rich or overly Ba-rich, which are HE 0305-4520 ([C/Fe] ~ 

+0.3, [Sr/Fe] 0.7, [Y/Fe] 0.4, [Ba/Fe] - +0.6) and 

HE 2156-3130 ([C/Fe] - +0.7, [Sr/Fe] 0.9, [Y/Fe] ~ 

—0.9, [Ba/Fe] ~ +0.5 ), the latter being the more metal-poor 
of the two. These stars perhaps warrant further study. 

We note that Figs. [22] and |24| show 6 stars with similar 
Sr/Ba ratios and similar C/H and Fe/H to the stars identified 
as s-process-rich stars. These stars also stand out clearly in 
a plot of [Ba/Fe] vs [C/Fe], see Fig. [28] occupying a similar 
region of the plot as the s-process-rich stars. As Eu has not 
been detected at the 3 sigma level according to our criteria they 
have not been classified in terms of the likely neutron-capture 
processes that have contributed to the heavy elements. It is in- 
teresting to look at these stars in more detail. We find that all 
these stars appear to be s-process-rich stars based on Eu abun- 
dances obtained, though we emphasise the abundances are not 
3 sigma detections and thus not completely reliable. HE 1430- 
1 123, the star with the lowest [Sr/Ba] of the six, and thus more 
closely associated with the s-II stars in Fig. [22] appears also 
to be a s-II star as we determine [Eu/Fe] ~ 1.4. The other five 
stars, HE 0231-4016, HE 0430-4404, HE 2150-0825, HE 2227- 
4044, and HE 2240-0412, which seem more closely associated 



Fig. 23. A plot of [Sr/Eu] vs [Fe/H] (upper panel) and [C/H] 
(lower panel). Symbols and lines have the same meanings as in 
Fig-E] with circles representing stars where Ba is undetected. 
The solar s-process value is 1.08. 
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Fig. 24. A plot of [Sr/Ba] vs [Ba/Fe]. Symbols and lines have 
the same meanings as in Fig. El w i m circles representing stars 
where Eu is undetected. 

with the remaining mildly Eu enhanced s-process-rich stars in 
Fig. 1221 appear to be similar stars with [Eu/Fe] ~ 0.7 <-> 1.0. 

Figure [28] also identifies a number of mildly carbon- 
enhanced stars without enhancement of neutron-capture ele- 
ments, a class of objects identified by Aoki et al. d2002ai. 
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Fig. 25. A plot of [Sr/Eu] vs [Eu/Fe]. Symbols and lines have Fig. 27. A plot of [Y/Ba] vs [Ba/Fe]. Symbols and lines have 
the same meanings as in Fig. [23] the same meanings as in Fig. 1261 



2 
1 

° n 
>- 

-1 

-2 



-4.0 -3.5 -3.0 -2.5 -2.0 
[Fe/H] 



-1 .5 



-1 .0 



2 
1 

-1 

-2 



-2 
[C/H] 



Fig. 26. A plot of [Y/Ba] vs [Fe/H] (upper panel) and [C/H] 
(lower panel). Symbols and lines have the same meanings as in 
Fig- El w i m circles representing stars where Eu is undetected. 



The five objects with [C/Fc] > 0.8 and [Ba/Fc] < 0.2 are, 
from largest [C/Fe] to smallest, HE 1351-1049, HE 1300-0641, 
HE 1330-0354, HE 1300-2201, and HE 1124-2335. 
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Fig. 28. A plot of [Ba/Fe] vs [C/Fe]. Symbols have the same 
meanings as in Fig. 1221 



6. Conclusions 

We have analysed snapshot spectra of a sample of 253 metal- 
poor stars, deriving abundances for 22 elements, where de- 
tectable, using an automated technique based on SME by 
Valenti & Piskunov ( 1996). The technique has been shown to 
give results in agreement with previous work within errors, 
when one considers differences in temperature scale. For our 
particular case (namely resolving power, spectral coverage and 
line list) the technique has been shown to be quite robust (at 
around the 0.1 dex level) even for S/N as low as 15. Our error 
estimates and comparisons with studies in the literature show 
the derived elemental abundances to be of moderate precision, 
relative and absolute errors of order 0.15 and 0.25 dex respec- 
tively. This work has dealt with the CH weak content of the 
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sample. The CH strong content will also be examined as part 
of the HERES programme. It would be interesting to obtain 
additional spectra for (the most interesting of) these stars, par- 
ticularly covering the Eu II lines at 6437 and 6645 A. 

Our main goal here has been to identify stars of interest, 
particularly those enhanced in r-process elements. We identi- 
fied 8 new r-II stars, 35 r-I stars and 3 s-process rich stars with 
strong Eu enhancement. Notably, the r-II stars were found only 
in a rather narrow metallicity range, —3.2 <j [Fe/H] ^ —2.6. 
These interesting stars should be investigated in more detail, 
and such work is underway as part of the HERES project large 
programme. The spectra obtained will be of much higher qual- 
ity, in terms of S/N, resolving power, and spectral coverage. 
The higher quality data, together with a careful manual spec- 
trum analysis, will naturally permit abundances and stellar pa- 
rameters to be obtained with higher precision and abundances 
for additional elements to be obtained. It will also be possible 
to obtain isotopic ratios for some elements. Ideally, such anal- 
yses should be performed with reference to deviations from the 
assumptions of LTE and ID model atmospheres. 

The results presented here provide a database of confirmed 
metal-poor stars, including a number of new r-process and s- 
process rich stars, which may be used for selection of stars 
for further studies. The results also provide a homogeneous 
database of moderate precision abundances for comparison 
with Galactic chemical evolution models. We stress that such 
comparisons must consider the limitations of our data set, most 
importantly the selection effects arising from the significant in- 
completeness of the data for some elements, but also the pre- 
cision of the abundances and the assumptions of the modelling 
such as LTE and the use of ID model atmospheres. 

We investigated trends and scatter in some measured abun- 
dance ratios of interest. Among the stars without strong C en- 
hancement, at about [Fe/H] J> —2.5 we find that the cos- 
mic scatter in any abundance ratio is small (perhaps even non- 
existent), implying that at around this level of enrichment the 
Galactic halo was reasonably well mixed. At lower metallic- 
ities C, Sr, Y, Ba and Eu, and perhaps Zr, show evidence for 
cosmic scatter, while the results for the Mg, Ca, Sc, Ti, Cr, Fe, 
Co and Ni still indicate at most small cosmic scatter within the 
errors of our analysis. Due to the difficulties in accurately es- 
timating errors, and thus disentangling the observational and 
modelling uncertainties from the real cosmic scatter, we con- 
clude that to determine the magnitude of the scatter in the cases 
where it is small will require high precision studies of a large 
number of stars. 
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Appendix A: Line List 

The lines and spectral windows used are listed in Tables |A~TT - 
IA.5I The most important atomic and molecular data, wave- 
length A, excitation potential x, and log gf, are tabulated. Line 
broadening data used may be obtained from the authors on re- 
quest. Below we briefly comment on line selection issues and 
data sources for each element. 

First, it is worth commenting on Si which is notably omit- 
ted from our analysis. Our spectra include two Si lines, at 
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3905 and 4102 A. The A4102 line lies in the wing of H<5 and 
thus presents a considerable challenge for automated analy- 
sis. The A3905 line is known to be blended by CH lines (e.g. 
Cayrel et al. 2004 1. Even with known CH lines from the lists 
of Plez & Jorrisen (see below) included in the synthesis, our Si 
abundances were found to differ from the literature values (see 
Sect. 13.5.21 . which are usually based on the A4102 line, and the 
difference was seen to correlate with carbon-enhancement. We 
have reasonable agreement with the literature for CS 22186- 
025 which has [C/Fe]^ —0.7, while we have mild disagree- 
ment for CS 31082-001 which has [C/Fe]~ 0.2 and strong dis- 
agreement for CS 22892-052 which has [C/Fe]^ 1. This sug- 
gests that the A3905 line may be affected by further blends of 
features involving C, and thus we chose to omit Si from our 
analysis. 

Carbon Abundances are derived from selected clean regions 
of CH A-X bands at 43 1 0^13 1 3 A (G band) and 4362-4367 A. 
We limited ourselves to these regions to reduce computing 
times due to the large number of components in these bands. 
Comparisons with results using larger regions showed no dif- 
ferences to those using the final adopted regions. We note, how- 
ever, that we consider it important to use at least two different 
regions of differing strengths, in particular because the G band 
region at 4310 A can become saturated in carbon-enhanced 
stars. The observations are not of sufficient quality to deter- 
mine isotopic abundances, and so we made the assumption that 
all carbon is in the form of 12 C and there is no 13 C. The adopted 
line list given in Tables |A~T1 and lA~2l was extracted from a list 
compiled by Plez & Jorissen (private communication), which 
is described by Hill et al. d2002b . 

Magnesium We employed 7 spectral features of Mg I. The 
features consist of the generally strong UV triplet near 3835 A 
from the 3 P°— 3 D transition, plus some weaker yet cleaner fea- 
tures farther to the red, the strongest of which is generally 
A4703. Our employed lines are practically the same as those 
adopted by Mc William et al. J1995bt and Norris et al. ( 1996) 
(see also Ryan et al. 1996), though we note that there are some 
significant differences in the adopted log gf values, particu- 
larly for the triplet where the two works adopt values which 
differ by as much as 0.5 dex. This warrants some discussion. 
Mc William et al. adopted log gf values for the triplet from 
the NIST compilation of Wiese et al. 1 1969), which are based 
on theoretical calculations of Weiss ( 1967) and re-normalised 
measurements by Penkin & Shabanova (1962); the compila- 
tion can be consulted for details. We note that the 3 P°— 3 D 
transition is actually a sextuplet, but due to the practically non- 
existent fine-structure splitting of the 3 D level only 3 lines are 
observed. Wiese et al. (1969) present log gf values for each 
of the 6 components, but McWilliam et al. seem only to have 
included the strongest component for each line. We have also 
adopted the NIST data, but have added together the appropriate 
components. Thus, two of the lines have log gf values of order 
0.1 dex stronger than those used by McWilliam et al dT995b I. 
Norris et al. seem not to have been aware of the NIST data, and 
adopted what amounts to empirical /-values guided by the de- 



sire to obtain abundances from these lines consistent with those 
from the other lines, particularly the A4703 line, where the 
theoretical calculations of Froese Fischer J1975i are employed 
(which we also adopt for these lines). Their adopted values for 
the triplet differ from our adopted values by as much as 0.6 dex. 
Thus, we would expect to find largely discrepant results from 
the triplet and the remaining lines with our /-values. However, 
in our sample we always found a good general global fit with a 
single abundance. We suspect this apparent discrepancy is due 
to the fact that we use a profile fitting technique with the line 
cores removed, while Norris et al. fit equivalent widths. The 
triplet lines are quite strong, having equivalent widths often of 
order 100 mA even for stars with [Fe/H]< —3, and our fits to 
the triplet show the observed cores to be typically much deeper 
than those from the best fit to the wings. This is a well known 
problem of LTE analyses with no modelling of chromospheres. 
This effect goes in the correct direction, in the sense that Norris 
et al. would need to underestimate the /-values in order to cor- 
rect for the overestimation of abundance due to the influence of 
the line cores. We note that, since the Norris et al. /-values for 
the triplet are calibrated to the weaker line abundances, their 
final abundances should be correct. 

Al U m i n i U m Only one line of Al I was employed, the resonance 
line A3961. The second resonance line A3944 was too blended 
to be considered reliable. Hyperfine structure was adopted from 
McWilliam et al. dl995bt . 

Calcium The Ca abundance is determined from 9 lines of 
Ca I. The /-values are adopted from the NIST critical com- 
pilations. 

Scandium The /-values are taken from Lawler & 
Dakin (1989). Hyperfine structure is taken directly from 
McWilliam et al. dl995bl ). 

Titanium Oscillator strengths for neutral lines are adopted 
from Blackwell et al. ( 1982a1 fT982bt . For singly ionised lines, 
the situation has been markedly improved by the experimental 
work of Pickering et al. J2001i . and their oscillator strengths 
have been adopted. For two useful lines where data were not 
available from this source, we adopted the solar values from 
Ryabchikova et al. ( 1994). We note that for the other Ti II lines 
we used these two sources agreed within 0. 1 dex. 

Vanadium The log gf data for V I is taken from Doerr 
et al. (1985), and for V II from Biemont et al. (1989). 
McWilliam et al. ( 1995b) noted the lack of hyperfine structure 
data for the transitions used here, and to our knowledge this sit- 
uation has not changed. However, as noted by McWilliam et al. 
these lines are weak and unsaturated and therefore our abun- 
dances should not be greatly affected by this omission. We note 
that in the case there was some effect, our abundances would 
be an upper limit. 
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Chromium The /-values are taken from the NIST compila- 
tion (Martin et al. 1988). The values are based on a critical 
averaging of number of sources, which are detailed in the com- 
pilation. 

Manganese We adopted oscillator strengths from Booth 
et al. ( 119841 ). The hyperfine structure was computed using the 
data compiled in Lefebvre et al. J2003> . which for the lines 
employed here draws on data from Davis et al. dl 97 1 I t. and 
Handrich et al. J1969i . as well as their own measurements. 

Iron A total of 55 lines were used of which 45 were Fe I lines 
and 10 Fe II. For the neutral lines the oscillator strengths of 
O' Brian et al. dl991> were adopted. Oscillator strengths for the 
Fe II lines were adopted from Schnabel et al. (2004) where 
possible. Data for some other useful lines were taken from 
Biemont et al. (119911 ) and Moity ( 119831 ). We note that the Fe II 
line /-values are of lower precision, typically 20%, than the 
Fe I lines, typically 7%. 

Cobalt Oscillator strengths for the employed lines were taken 
from Nitz et al. ( 1999) where possible, otherwise from Cardon 
et al. J1982i . Hyperfine structure was included using the data 
of Pickering dl996i . 

Nickel Four lines of Ni I are employed. For 3 lines we em- 
ployed the oscillator strengths from Blackwell et al. d!989i ; 
for the remaining line we used the /-value from Huber & 
Sandeman (fT980l . 

Zinc Two neutral lines were employed and the /-values taken 
from Biemont & Godefroid fTWfl . 

Strontium One line of Sr I and two of Sr II were analysed. 
Hyperfine and isotopic splitting was accounted for in the A4215 
line, using the data from Borghs et al. ( 1983 ), where solar iso- 
topic ratios have been assumed. There are unfortunately no data 
to our knowledge for the upper state of the other resonance 
line A4077. However, noting that the splitting in the A4215 
line is dominated by the common lower level 2 <Si/2, we com- 
puted hyperfine splitting for the A4077 line where we neglected 
the splitting of the upper level. The main isotope shifts for the 
A4215 line are interpreted as due to the volume effect which 
affects the s state most strongly (e.g. Cowan [T981I pg 506). 
We therefore adopted the same isotopic shifts for A4077 as for 
the A4215 line. We did not find hyperfine structure data for the 
remaining line. However, as it is weak, this should not be im- 
portant. 

Yttrium The /-values for the employed lines of Y II are from 
Hannaford et al. ( 1982). 89 Y has a small nuclear spin of / = 
1/2, and thus hyperfine splitting is negligible. 



Zirconium Three lines of Zr II are analysed employing the /- 
values from Biemont et al. ( 1 98U . 

Barium Our analysis of Ba is based on 2 lines of Ba II, with 
oscillator strengths from Gallagher ( 119671 ) which for common 
lines are in excellent agreement with Davidson et al. {T992). 
Hyperfine structure was included from Mc William ( 1998 ), as- 
suming the pure r-process isotopic composition from that pa- 
per. As pointed out by Mc William ( 1998 ), if a star would have 
a significant s-process contribution, then our derived Ba abun- 
dance would be a lower limit. This is because the s-process 
causes a higher fraction of the Ba to be in even isotopes which 
have no hyperfine splitting, and thus the stronger lines saturate 
more quickly, meaning the same line strength would require a 
higher abundance. 

Lanthanum The oscillator strengths of Lawler et al. (120011 
are used. Hyperfine structure is included using the constants 
from Lawler et al., which for our lines draws from their own 
measurements and those of Hohle et al. (I1982> . 

Cerium We employed 10 lines of Ce II, the same as those 
employed by Sneden et al. ( 1996 ), with the addition of the line 
A3999. We adopted the same oscillator strengths as used by 
Sneden et al (1996). For A3999 we adopted the VALD value, 
which is essentially a solar value averaged from various sources 
(see Piskunov et al.[l995 ). 

Neodymium We employ 9 lines of Nd II using laboratory /- 
values from Den Hartog et al. (j2003 >, which are a significant 
improvement on those previously available. 

Samarium Eight lines of Sm II are employed and the /- 
values from Biemont et al. (Q989) were used. 

Europium Our Eu abundances are based on 4 lines of Eu II. 
Lines at 4435 and 4522 A were found to be too strongly 
blended in non-r-process-enhanced stars to be employed. The 
line A3971 was considered too blended for automated anal- 
ysis. Oscillator strengths and hyperfine splitting from Lawler 
et al. (2001 ) were employed, assuming solar r-process isotopic 
fractions (supported by measurements by Sneden et al. 12002b . 

Appendix B: Error Estimates 

Our method for estimating error propagation is similar to the 
approach devised by Mc William et al. ( 1995b), but modified 
to suit our abundance analysis method. The important differ- 
ence is that our approach fits the spectral features of a given 
element globally rather than fitting individual lines. Further, as 
discussed in Sect. 13.41 we will make a distinction between ab- 
solute (the uncertainty in the absolute abundance) and relative 
(the uncertainty in the relative abundance between stars) error 
estimates. First, we develop the formalism quite generally, and 
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Table A.l. Molecular line data for A-X system of the CH 
molecule near 4310 A from Plez & Jorissen. The spectral win- 
dows employed for this molecular band are defined with re- 
spect to an arbitrary A and span from A + dAbiue to A + dA re d. 



Species 


A 


X 


logs/ 




[A] 


[eV] 




CH 


4310.150 






CH 


4310.038 


0.096 


-3.183 


CH 


4310.090 


0.096 


-1.412 


CH 


4310.110 


0.096 


-1.474 


CH 


4310.149 


0.431 


-3.202 


CH 


4310.162 


0.096 


-3.025 


CH 


4310.203 


0.431 


-1.443 


CH 


4310.220 


0.431 


-1.505 


CH 


4310.272 


0.431 


-3.039 


CH 


4310.404 


0.096 


-3.183 


CH 


4310.430 


0.408 


-1.581 


CH 


4310.458 


0.096 


-1.412 


CH 


4310.504 


0.432 


-3.202 


CH 


4310.508 


0.408 


-2.966 


CH 


4310.556 


0.432 


-1.443 


CH 


4310.679 


0.072 


-1.550 


CH 


4310.729 


0.408 


-1.581 


CH 


4310.757 


0.072 


-2.952 


CH 


4310.807 


0.408 


-2.966 


CH 


4310.811 


0.408 


-3.155 


CH 


4310.889 


0.408 


-1.508 


CH 


4310.937 


0.388 


-1.677 


CH 


4310.969 


1.107 


-2.269 


CH 


4310.991 


0.072 


-1.550 


CH 


4311.045 


0.388 


-2.879 


CH 


4311.069 


0.072 


-2.952 


CH 


4311.075 


0.408 


-3.155 


CH 


4311.075 


0.072 


-3.136 


CH 


4311.084 


0.987 


-4.123 


CH 


4311.145 


0.987 


-4.123 


CH 


4311.153 


0.408 


-1.508 


CH 


4311.153 


0.388 


-1.677 


CH 


4311.153 


0.072 


-1.477 


CH 


4311.196 


1.108 


-3.922 


CH 


4311.196 


1.108 


-2.229 


CH 


4311.261 


0.388 


-2.879 


CH 


4311.328 


0.371 


-1.805 


CH 


4311.348 


0.072 


-3.136 


CH 


4311.394 


0.388 


-3.105 


CH 


4311.426 


0.072 


-1.477 


CH 


4311.473 


0.371 


-1.805 


CH 


4311.476 


0.371 


-2.771 


CH 


4311.502 


0.388 


-1.587 


CH 


4311.502 


0.051 


-1.646 


CH 


4311.545 


1.108 


-2.269 


CH 


4311.547 


0.347 


-2.386 


CH 


4311.567 


0.357 


-2.002 


CH 


4311.580 


0.388 


-3.105 


CH 


4311.592 


0.347 


-2.386 


CH 


4311.612 


0.051 


-2.865 


CH 


4311.618 


0.662 


-2.054 


CH 


4311.618 


0.662 


-2.054 


CH 


4311.621 


0.371 


-2.771 


CH 


4311.655 


0.357 


-2.002 


CH 


4311.688 


0.388 


-1.587 


CH 


4311.722 


1.108 


-3.923 


CH 


4311.727 


0.051 


-1.646 


CH 


4311.729 


1.108 


-2.229 


CH 


4311.776 


0.357 


-2.632 


CH 


4311.837 


0.051 


-2.866 


CH 


4311.861 


0.347 


-2.432 


CH 


4311.863 


0.357 


-2.632 


CH 


4311.897 


0.371 


-3.052 


CH 


4311.906 


0.347 


-2.432 


CH 


4311.978 


0.051 


-3.085 


CH 


4312.017 


0.371 


-3.052 


CH 


4312.045 


0.371 


-1.688 


CH 


4312.088 


0.051 


-1.557 


CH 


4312.153 


0.033 


-1.773 


CH 


4312.164 


0.371 


-1.688 


CH 


4312.172 


0.051 


-3.085 


CH 


4312.280 


0.051 


-1.557 


CH 


4312.304 


0.033 


— 1.773 


CH 


4312.304 


0.033 


-2.759 


CH 


4312.317 


0.358 


-3.0O7 


CH 


4312.386 


0.358 


-3.0O7 


CH 


4312.456 


0.033 


-2.759 


CH 


4312.527 


0.358 


-1.831 


CH 


4312.594 


0.358 


-1.831 


CH 


4312.594 


0.019 


-1.969 


CH 


4313.620 






CH 


4313.377 


0.020 


-2.985 


CH 


4313.446 


0.020 


-2.985 


CH 


4313.590 


0.020 


-1.800 


CH 


4313.660 


0.020 


-1.800 


CH 


4313.876 


0.009 


-3.033 


CH 


4313.906 


0.009 


-3.033 



dA bluo 

[mA] 



<lA rcd 
[mA] 



Table A.2. Molecular line data for A-X system of the CH 
molecule near 4363 A from Plez & Jorissen. The spectral win- 
dow employed for this molecular band is defined with respect 
to two arbitrary wavelengths A and span from A + dAbiue to 
A + rfA rc d- 



Species 


\ 


X 


l°g £// 




rXi 


[eV] 




CH 


4363. jUU 






CH 


43oZ.UZl 


1 .247 


— 4.054 


CH 


A 1 £0 f\H 1 

43oZ.Uol 


1.247 


— 1 .976 


CH 


43oZ. 1 12. 


1 .247 


— 1.943 


CH 


43oZ.ZUZ 


0.777 


— 1 .982 


CH 


43oZ.Z()4 


1 .247 


— 4. 149 


CH 


43oZ.Zdd 


0.777 


— 3.284 


CH 


4ioZ.3il 


0.777 


— 1.917 


CH 


43oZ.54y 


0.557 


— 1.749 


CH 


4362.551 


0.557 


— 3.392 


CH 


43oZ.oy / 


0.777 


— 1 .982 


CH 


Ali^ n At\ 

4 JoZ. /4y 


0.557 


— 1.705 


CH 


4ioZ. /dU 


0.777 


— 3.285 


CH 


43oZ.y&3 


0.777 


— 1.917 


CH 


43o3.Uno 


0.227 


— 3.380 


CH 


43o3.U8 / 


0.227 


— 1.716 


CH 


43o3.1oz 


1 .248 


— 4.055 


CH 


43o3. 18o 


1 .248 


— 1 .976 


CH 


43o3.Z4Z 


1 .248 


— 1.943 


CH 


43o3.Z/o 


1 .248 


— 4. 149 


CH 


43o3. Zy 1 


0.227 


— 1.673 


CH 


4363.461 


0.558 


— 1.749 


CH 


4363.463 


0.558 


— 3.392 


CH 


43oi.oU3 


0.558 


— 1 .705 


CH 


Alii A f\1 A 

43o4.U34 


0.228 


— 3.380 


CH 


43o4.U3o 


0.228 


— 1.716 


CH 


Alii A l r A 

43o4. 1 14 


1 .289 


— 4.093 


CH 


4364. 167 


1 .289 


— 1 .947 


CH 


Alii A 10 1 

43o4.ini 


0.228 


— 1.673 


CH 


4364.253 


1 .289 


— 1.916 


CH 


43o4.Zy3 


1 .289 


— 4. 1 8 1 


CH 


4365.416 


1 .290 


— 4.093 


CH 


43o3.43U 


1 .290 


— 1 .947 


CH 


A^az Ann 
430D.4/ / 


1 .290 


— 1.916 


CH 


43o3.3ZZ 


1 .290 


— 4. 1 8 1 


CH 


43oD.yzy 


0.987 


— 3.729 


CH 


4300.U1U 


0.987 


— 3.729 


CH 


43oo.Z3U 


0.987 


— 3.828 


CH 


4366.312 


0.987 


-3.828 


CH 


4366.398 


0.597 


-3.426 


CH 


4366.407 


0.597 


-1.699 


CH 


4366.482 


0.268 


-3.414 


CH 


4366.497 


0.268 


-1.666 


CH 


4366.520 


1.333 


-4.129 


CH 


4366.573 


0.597 


-1.659 


CH 


4366.588 


1.333 


-1.920 


CH 


4366.647 


1.333 


-1.892 


CH 


4366.663 


0.268 


- 1 .627 


CH 


4366.699 


1.333 


-4.211 


CH 


4366.897 


0.805 


- 1 .903 


CH 


4366.928 


0.987 


-3.495 


CH 


4366.930 


0.805 


-3.326 


CH 


4366.991 


0.987 


-3.495 


CH 


4367.172 


0.805 


-1.848 



d^bluo 
[mA] 



[mA] 



later specify the difference between the calculations for these 
two types of errors. 

We consider the propagation of errors in model atmosphere 
parameters T D g , log g, £. For low metallicity models of the type 
used here, typical errors in metallicity have negligible effect on 
the model structure and can be neglected. Errors due to prop- 
agation of uncertainties in logg/, observational error, contin- 
uum placement and spectrum modelling uncertainties are also 
considered. For compactness, in the following discussion we 
define the abundance parameter e = log e. Assuming the be- 
haviour of the considered abundance with small changes in 
these parameters can be approximated by a first-order Taylor 
expansion, i.e. that the abundance e = loge varies approxi- 
mately linearly with changes in parameters on scales of the pa- 
rameter uncertainties, we obtain (e.g. Mc William et al. 1995b 
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Table A.3. Data for atomic lines for species lighter than Fe. 
The spectral window employed for the line is from A + dAbiue 
to A+dA re( j- Lines which are removed in carbon-enhanced stars 
are marked with asterisks in the wavelength column. 



Species 


A 


X 


logs/ 






Refs. 




[A] 


[eVj 




[mA] 


[mA] 




Mgl 


3829.355 


2.707 


-0.208 


-150 


100 


WSM69 


Ma I 


3832.304 


2.710 


0.270 


-300 


200 


WSM69 


Mg 1 


3838.292 


2.715 


0.490 


-230 


150 


WSM69 


Mgl 


•4057.505 


4.346 


-0.890 


-120 


100 


FF75 


Me 1 


•4167 27 1 


4 346 


Q 7 [Q 


250 


mo 


FF75 


Mg 1 


4571.096 


0.000 


-5.393 


-300 


300 


WSM69 


Mgl 


4702.990 


4.330 


-0.380 


-300 


300 


FF75 


Al I 


3961.529 


0.014 


-0.336 


-200 


170 


WM80 


Ca 1 


4226 728 


000 


244 


— 150 


400 


WM80 


Ca 1 


•4283.011 


1.886 


-0.224 


-120 


300 


WM80 


Ca 1 


•4289.367 


1.879 


-0.304 


-120 


170 


WM80 


Ca 1 


4318.652 


1.890 


-0.207 


-120 


120 


WM80 


Ca 1 


•4425.437 


1.879 


-0.358 


-300 


300 


WM80 


Ca 1 


•4434.957 


1.886 


-0.005 


-300 


50 


WM80 


Ca 1 


4454.780 


1.898 


0.258 


-50 


200 


WM80 


Ca 1 


4455.887 


1.899 


-0.526 


-300 


300 


WM80 


Ca 1 


4578.551 


2.521 


-0.558 


-200 


loo 


WM80 


Sc 11 


4246.837 


0.315 


0.240 


-350 


250 


LD89 


Sc 11 


•4314.095 


0.618 


-0.100 


-250 


70 


LD89 


Sc 11 


•4400.399 


0.605 


-0.540 


-150 


200 


LD89 


Sc 11 


4415.563 


0.595 


-0.670 


-180 


300 


LD89 


Sc 11 


4670.417 


1.357 


-0.580 


-150 


50 


LD89 


Ti 1 


3989.759 


0.021 


-0.140 


-200 


100 


GBP89, BPSL82 


n i 


3998.636 


0.048 


0.000 


-250 


150 


GBP89, BPSL82 


Ti I 


4533 749 


848 


530 


— 150 


300 


GBP89 BMPS82 


Ti 1 


4534.776 


0.836 


0.340 


-300 


300 


GBP89, BMPS82 


Ti 1 


4656.469 


0.000 


-1.290 


-300 


300 


GBP89, BPSL82 


Ti 1 


•4681.909 


0.048 


-1.020 


-300 


150 


GBP89, BPSL82 


Ti 11 


•4337 915 


1 (ISO 


—0 960 


— 150 


1 *i0 


PTP01 


Ti 11 


•4394.051 


1.221 


-1.780 


-100 


300 


PTP01 


Ti 11 


•4395.850 


1.243 


-1.930 


-200 


100 


PTP01 


Ti 11 


•4417.719 


1.165 


-1.190 


-200 


300 


PTP01 


Ti 11 


4443.794 


1.080 


—0.720 


— 300 


250 


PTP01 


Ti 11 


4468.507 


1.131 


-0.600 


-150 


300 


RHLPS94 


Ti 11 


•4470.857 


1.165 


-2.020 


-250 


150 


PTP01 


Ti 11 


4501.273 


1.116 


-0.770 


-300 


300 


PTP01 


Ti 11 


4563.761 


1.221 


-0.690 


-300 


300 


PTP01 


Ti 11 


4571.968 


1.572 


-0.320 


-300 


150 


PTP01 


Ti 11 


4589.958 


1.237 


— 1.620 


—300 


300 


RHLPS94 


V 1 


•4379.230 


0.301 


0.550 


-150 


150 


DKKWZ85 


V 1 


•4389.976 


0.275 


0.270 


-150 


150 


DKKWZ85 


V 11 


•3951.960 


1.480 


-0.784 


-150 


60 


BGFML89 


Cl 1 


4254.332 


0.000 


-0.114 


-150 


200 


see MFW88 


Cl 1 


•4274.796 


0.000 


-0.231 


-200 


100 


see MFW88 


Cr 1 


•4289.716 


0.000 


-0.361 


-150 


100 


see MFW88 


Mil 1 


4030.763 


0.000 


-0.470 


-150 


250 


BBPS84 


Mil 1 


4033.060 


0.000 


-0.618 


-200 


150 


BBPS84 


Mil 1 


•4034.492 


0.000 


-0.811 


-150 


300 


BBPS84 


Mil 1 


4754.040 


2.282 


-0.086 


-200 


200 


BBPS84 


Mnl 


4823.496 


2.319 


0.144 


-100 


200 


BBPS84 



References: BBPS84 = Booth et al. 1984 . BGFML89 = Biemont et al. 1989 . BMPS82 = Blackwell 
et al. I1982al. BPSL82 = Blackwell et al . I1982bl . DKKWZ85 = Doerr et al. [79851 , FF75 = Froese 
Fischer 119751. GBP 89 = Grevesse et al. 119891. LD 89 = Lawler & Dakin 119891 , M FW88 = Martin, 
Fuhr & Wiese 119881. P TP01 = Pickering etalJMOl , RHLPS94 = Ryabchikova et al. [79941 , WM80 = 
Wiese & Martin 1980 , WSM = Wiese et al. 119691. 



Taylor 19821 for the variance in the abundance 

de 



de_ 
df 

2 



-a 
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de 



.dloggr, 

of (obs) + of (cont) + of (model) 
de de de 



dTdlogg 
de de 
'Moi^9| CTlos9 ' ? 



dT d£ 



(B.l) 



where of is the variance in parameter i, and oy the covari- 
ance of i and j. The variances of (obs) and of (cont) represent 
the variance in the abundance e due to observational error and 



Table A.4. Data for atomic lines of Fe. The spectral window 
employed for the line is from A + dAbiuc to A + <iA re d- Lines 
which are removed in carbon-enhanced stars are marked with 
asterisks in the wavelength column. 



Species 


X 


X 


logs/ 


dA bluo 


dA„ d 


Refs. 




[A] 


[eVJ 




[mA] 


[mA] 




Fel 


•3856.372 


0.052 


-1.286 


-200 


400 


OWLWB91 


Fe 1 


•3859.911 


0.000 


-0.710 


-400 


350 


OWLWB91 


Fe 1 


3865.523 


1.011 


-0.950 


-250 


250 


OWLWB91 


Fe 1 


•3878.018 


0.958 


-0.896 


-300 


100 


OWLWB91 


Fe 1 


4005.242 


1.557 


-0.583 


-150 


150 


OWLWB91 


Fe 1 


4045.812 


1.485 


0.284 


-100 


150 


OWLWB91 


Fe 1 


4063.594 


1.558 


0.062 


-100 


300 


OWLWB91 


Fe 1 


4071.738 


1.608 


-0.008 


-100 


300 


OWLWB91 


Fe 1 


4114.445 


2.831 


-1.303 


-250 


250 


OWLWB91 


Fe 1 


4132.058 


1.608 


-0.675 


-300 


100 


OWLWB91 


Fe 1 


4143.868 


1.557 


-0.511 


-200 


300 


OWLWB91 


Fe 1 


•4175.636 


2.845 


-0.827 


-150 


300 


OWLWB91 


Fe 1 


•4187.039 


2.449 


-0.514 


-150 


200 


OWLWB91 


Fe 1 


•4199.095 


3.047 


0.156 


-200 


300 


OWLWB91 


Fe 1 


4202.029 


1.485 


-0.689 


-200 


150 


OWLWB91 


Fe 1 


4222.213 


2.449 


-0.914 


-300 


200 


OWLWB91 


Fe 1 


4227.426 


3.332 


0.266 


-150 


150 


OWLWB91 


Fe 1 


•4233.603 


2.482 


-0.579 


-150 


100 


OWLWB91 


Fe 1 


•4250.120 


2.469 


-0.380 


-100 


300 


OWLWB91 


Fe 1 


4250.787 


1.557 


-0.713 


-300 


300 


OWLWB91 


Fe 1 


4260.474 


2.400 


0.077 


-150 


200 


OWLWB91 


Fe 1 


•4271.154 


2.449 


-0.337 


-300 


100 


OWLWB91 


Fe 1 


•4271.761 


1.485 


-0.173 


-150 


300 


OWLWB91 


Fe 1 


4282.403 


2.176 


-0.779 


-250 


150 


OWLWB91 


Fe 1 


•4325.762 


1.608 


0.006 


-250 


100 


OWLWB91 


Fe 1 


•4375.930 


0.000 


-3.031 


-150 


300 


OWLWB91 


Fe 1 


•4383.545 


1.485 


0.208 


-300 


300 


OWLWB91 


Fe 1 


•4404.750 


1.557 


-0.147 


-250 


200 


OWLWB91 


Fe 1 


•4415.123 


1.608 


-0.621 


-300 


200 


OWLWB91 


Fe 1 


•4430.614 


2.223 


-1.728 


-250 


250 


OWLWB91 


Fe 1 


•4442.339 


2.198 


-1.228 


-300 


300 


OWLWB91 


Fe 1 


4443.194 


2.858 


-1.043 


-100 


250 


OWLWB91 


Fe 1 


4447.717 


2.223 


-1.339 


-300 


300 


OWLWB91 


Fe 1 


4489.739 


0.121 


-3.899 


-300 


200 


OWLWB91 


Fe 1 


4494.563 


2.198 


-1.143 


-300 


300 


OWLWB91 


Fe 1 


4528.614 


2.176 


-0.887 


-100 


300 


OWLWB91 


Fe 1 


4602.941 


1.485 


-2.208 


-300 


300 


OWLWB91 


Fe 1 


•4736.773 


3.211 


-0.752 


-200 


200 


OWLWB91 


Fe 1 


4872.137 


2.882 


-0.567 


-300 


300 


OWLWB91 


Fe 1 


4890.755 


2.876 


-0.394 


-350 


350 


OWLWB91 


Fe 1 


4891.492 


2.852 


-0.112 


-300 


300 


OWLWB91 


Fe 1 


4918.994 


2.845 


-0.342 


-350 


350 


OWLWB91 


Fe 1 


4920.503 


2.832 


0.068 


-400 


100 


OWLWB91 


Fe 1 


4938.814 


2.875 


-1.077 


-300 


200 


OWLWB91 


Fe 1 


4939.687 


0.859 


-3.252 


-250 


300 


OWLWB91 


Fe 11 


•4178.862 


2.583 


-2.443 


-100 


200 


SSK04 


Fe 11 


4233.172 


2.583 


-1.809 


-150 


200 


SSK04 


Fe 11 


4416.828 


2.778 


-2.540 


-240 


300 


M83 


Fe 11 


4508.289 


2.856 


-2.318 


-300 


300 


BBKAP9I 


Fe 11 


4515.343 


2.844 


-2.362 


-100 


300 


SSK04 


Fe 11 


4520.224 


2.807 


-2.550 


-200 


300 


M83 


Fe 11 


•4541.524 


2.856 


-2.990 


-100 


200 


M83 


Fe 11 


•4555.893 


2.828 


-2.250 


-150 


100 


SSK04 


Fe 11 


4583.839 


2.807 


-1.740 


-250 


300 


SSK04 


Fell 


4923.927 


2.891 


-1.206 


-400 


400 


SSK04 



References: BBKAP91 = Biemo nt et al. 179911 . M83 = Moity [79831 . OWLWB91 = O'Brian et al JiWfl 
SSK04 = Schnabel et al. 120041 



continuum placement uncertainties respectively. The variance 
of (model) represents the variance in the abundance e due to 
spectrum modelling uncertainties, such as the assumptions of 
ID modelling and LTE. Terms involving covariances of inde- 
pendent parameters have been omitted. We have used the fact 
that de/d\og gf = 1. 

We adopt ax,£ = and a\ ogg ^ = —0.02 following 
McWilliam et al. (1995bi. We obtained similar estimates for 
these quantities from our own numerical experiments for a sub- 
sample of stars representative of the complete sample. Based 



on these experiments we adopted OT,iogg 
expression for the total variance becomes 



-22. Thus, the 



de\ 2 f de 

of) aT+ UW aioS9 

+^lo Sg f + <?e(obs) + of (cont) 



' de s 
of (model) 



(B.2) 
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J de de de de 1 

+2 \dfd\ogg (JT ' log9 + d^Y^^j ■ 

The required partial derivatives are obtained individually for 
each star and element by direct determination of abundances 
for shifted stellar parameters (see Sect. 13. 3> . 

Now it simply remains to specify each variance (or stan- 
dard deviation). The input variances, however, will depend on 
whether we wish to estimate the absolute error or the relative 
error. 



the error due to uncertainty in the continuum placement can be 
approximated by cr £ (cont) s» cr e (obs) j\fm, where the 
factor accounts for the fact that we typically have m times more 
independent pixels to define the continuum than the line. Based 
on inspection of spectra we adopt m = 5 for all elements ex- 
cept carbon, where we adopt m = 1 since it is determined from 
wide molecular bands. The estimates for m approximately ac- 
count for possible correlations between pixels introduced by 
very weak lines. The uncertainty in log g/ will not contribute 
to relative errors, and thus we adopt <J\ ogg f = 0. 



B. 1. Relative error estimates 



B.2. Absolute error estimates 



The relative error estimates are of interest for comparison of 
abundances within the sample. In particular, we wish to es- 
timate the amount of abundance scatter which should be at- 
tributed to uncertainties in the data and analysis. 

First, from Sect. 12. II we adopt ot = 100 K. If the measure- 
ment errors in the observed data are normally distributed, the 
formal random errors in derived spectrum modelling parame- 
ters due to this measurement error are given by the SME opti- 
misation routine from the estimated covariance matrix, which 
is the inverse of the curvature matrix at the solution (see e.g. 
Press et al l 1992b . These errors represent the propagation of the 
random observational uncertainties in each point in the spec- 
trum (dominated by photon noise) to the relevant parameter and 
is computed by the parameter optimisation procedure. Such 
contributions are obtained for cr e (obs), <J\ ogg and a^; these 
random components of a\ ogg and are usually quite small, 
around 0.05 dex and 0.05 km s _1 , while a e (obs) varies signif- 
icantly from element to element. To this random component of 
Clog g we add 0.22 dex reflecting the error due to uncertainty in 
T e ff , and to the random component of <jj we add 0. 1 km s _1 re- 
flecting the error due to uncertainty in log g. These numbers are 
based on the numerical experiments for a representative sub- 
sample which were mentioned above. 

As mentioned above, errors arise due to uncertainties in the 
spectrum modelling, for example, assumptions of ID geometry 
or LTE. The contributions of these errors to relative errors in 
abundances would be expected to cancel if all the stars were 
identical; however, across our sample, which covers a wide 
range of stellar atmosphere parameters, we can only expect par- 
tial cancellation. Thus, such uncertainties will lead to a degree 
of scatter in any quantity derived from the spectrum, arising 
from differences in modelling errors from star to star. The rel- 
ative component of a e (model) has been estimated at 0.1 dex. 
This is simply an order-of-magnitude estimate based on indica- 
tions from 3-D modelling, (e.g. Asplund 2004 2005 ) and non- 
LTE spectrum modelling (e.g. Korn et al. 2003 Asplund 2005 1. 
However it should be noted this quantity should vary with el- 
ement and employed spectral feature. We have neglected any 
relative error in log g and £ which might arise from modelling 
uncertainties. 

Since the relevant continuum points are located close to 
lines employed in the analysis, they will have very similar er- 
ror bars to the points in the line, and thus the sensitivity of the 
abundance to continuum placement errors will be essentially 
identical to that for the points in the line. Thus we expect that 



While relative errors are generally of most interest, it is also 
important to have some estimate of the absolute error in our 
obtained abundances. The calculation of absolute errors fol- 
lows that of the relative errors with a few changes which we 
now list. 

The main differences arise in the modelling uncertainties, 
as there is no cancellation as in the relative error case. To 
oiogg and cj£ we add an additional 0.1 dex and 0.1 km s _1 
respectively to account for modelling uncertainties. We esti- 
mate er e (model) at 0.15 dex. Further, the errors cr\ ogg f must 
be included in the absolute errors. Since our procedure glob- 
ally fits the spectral lines of a given element, the value for 
Clog g f should be representative of the typical error for the cho- 
sen spectral lines. In Table lB.il we provide an estimated aver- 
age uncertainty for each element with reference to the original 
literature, which are adopted for &\ ogg f. 

B.3. Errors in abundance ratios 

We are often interested in elemental abundance ratios. 
Following Mc William et al. ( 1995b I, if the abundances ea and 
Eb are expressed in logarithms such that £a/b = £ A — £b, the 
variance in the abundance ratio £a/b 1S given by 



°A/B = G 



<J B — 2(74 



(B.3) 



where a\ and a B are the variance in ea and eb respectively 
and ga.b is the covariance of ea and eb given by 



d£ A 9e b 
dT dT 



3ea 9e 
d log g d log g 



B ^ 2 

a log3 



Oea ds 



B ^ 2 



ds A Qe b 
dT dlogg 
( ds A d£ B 



+ 



\ <9£ dlogg + dlogg <9£ 



Pea d£ B 
dlogg dT 

8e a 9e b 



or. log 



Clog g,£- 



(B.4) 



Thus, as pointed out by Mc William et al. ( 1995b I, there may 
be partial cancellation of errors if the abundances of elements 
A and B have similar sensitivity to atmospheric parameters, 
or partial compounding of errors if the elements have contrary 
sensitivity. Note, abundance ratio error estimates may be com- 
puted in both the absolute and relative senses. 
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Table A.5. Data for atomic lines of species heavier than Fe. 
The spectral window employed for the line is from A + dAbiue 
to A+dA rc d- Lines which are removed in carbon-enhanced stars 
are marked with asterisks in the wavelength column. 



Species 


A 


X 


log 9/ 






Refs. 




rAl 

t rt l 


Lev j 




ImAl 






Col 


•3842.046 


0.923 


-0.770 


-200 


100 


CSSTW82 


Co I 


•3845461 


0.923 


0.010 


-150 


1511 


CSSTW82 


Co I 


3894.073 


1.049 


0.090 


-50 


200 


NKWL99 


Co I 


3995 30^ 


1 


140 


250 






Co I 


4118.767 


1.049 


-0.470 


-100 


200 


NKWL99 


Co I 


4121.311 


0.922 


-0.300 


-250 


250 


NKWL99 


Ni I 


3775.565 


0.423 


-1.408 


-200 


200 


BBPL89 


Nil 


3783.524 


0.423 


-1.304 


-70 


300 


HS80 


Ni I 


3807.138 


0.423 


-1.220 


-200 


150 


BBPL89 


Ni I 


3858 797 


4"* 1 


951 


300 


1 ( in 


RRPI RQ 


Zn 1 


•4722.163 


4.030 


-0.390 


-300 


300 


BG80 


Zn 1 


4810.537 


4.078 


-0.170 


-300 


300 


BG80 


St 1 


•4607.327 


0.000 


0.283 


-300 


100 


MB87 


SrD 


4077.709 


0.000 


0.158 


-150 


100 


PBL95 


Sr 11 


4215 519 


000 


q I55 


300 




PBL95 


Y 11 


3774.331 


0.130 


0.210 


-250 


150 


HLGBW82 


Y 11 


•3788.694 


0.104 


-0.070 


-150 


300 


HLGBW82 


Y 11 


3818.314 


0.130 


-0.980 


-300 


150 


HLGBW82 


Y 11 


3950.352 


0.104 


-0.490 


-200 


300 


HLGBW82 


Y 11 


•4398.013 


0.130 


-1.000 


-250 


150 


HLGBW82 


Y 11 


4883 684 


1 084 


070 


200 




1 II .1 It) W oZ 


ZrH 


•4161.213 


0.713 


-0.720 


-300 


130 


BGHL81 


Zrll 


4208.985 


0.713 


-0.460 


-200 


200 


BGHL8 1 


ZrH 


4317.299 


0.713 


-1.380 


-300 


300 


BGHL81 


Ba 11 


4130.645 


2.722 


0.560 


-150 


300 


G67 


Ba 11 


4554.000 


0.000 


0.163 


-500 


220 


G67 


La 11 


•3988.515 


0.403 


0.210 


-300 


300 


LBS0I 


La 11 


3995.745 


0.173 


-0.060 


-250 


150 


LBS0I 


La 11 


40R6 70Q 


000 


070 


200 


mo 


LBS01 


La 11 


4123.218 


0.321 


0.130 


-300 


180 


LBS0I 


La 11 


•4322.503 


0.173 


-0.930 


-300 


150 


LBS0I 


La 11 


•4333.753 


0.173 


-0.060 


-300 


100 


LBS0I 


Cell 


3999.237 


0.295 


0.232 


-100 


300 


VALD95 


Co 11 


4073.474 


0.478 


0.320 


-150 


150 


SMPCBA96 


Ce 11 


•4083.222 


0.701 


0.240 


-150 


200 


SMPCBA96 


Cell 


4120.827 


0.320 


-0.240 


-300 


300 


SMPCBA96 


Cell 


•4127.364 


0.684 


0.240 


-300 


70 


SMPCBA96 


Ce 11 


4222 597 


172 


180 


j 50 


1 ^0 


GS94 


Ce 11 


•4418.780 


0.864 


0.310 


-200 


150 


SMPCBA96 


Cell 


4486.909 


0.295 


-0.360 


-300 


300 


GS94 


Cell 


4562.359 


0.478 


0.330 


-300 


300 


GS94 


Ce 11 


4628.161 


0.516 


0.260 


-300 


300 


GS94 


Ndll 


•4018.823 


0.064 


-0.850 


-300 


120 


DLSC03 


Ndll 


4021.327 


0.321 


-0.100 


-200 


200 


DLSC03 


Ndll 


4061.085 


0.471 


0.550 


-300 


300 


DLSC03 


Nd II 


■4V JU V._U J 


1) 064 


570 


200 


illll 


ni srn3 


Ndll 


•4109.448 


0.321 


0.350 


-150 


150 


DLSC03 


Ndll 


4232.374 


0.064 


-0.470 


-200 


200 


DLSC03 


Ndll 


•4358.161 


0.321 


-0.160 


-200 


150 


DLSC03 


Ndll 


4446.384 


0.205 


-0.350 


-300 


300 


DLSC03 


Ndll 


•4462.979 


0.559 


0.040 


-300 


300 


DLSC03 


Smll 


3896.972 


0.041 


-0.578 


-100 


200 


BGHL89 


Sm II 


4068.324 


0.434 


-0.710 


-150 


150 


BGHL89 


Smll 


•4318.936 


0.277 


-0.270 


-150 


150 


BGHL89 


Sm 11 


4499.475 


0.248 


-1.010 


-200 


200 


BGHL89 


Smll 


4519.630 


0.544 


-0.432 


-300 


300 


BGHL89 


Sm 11 


4537.954 


0.485 


-0.230 


-300 


300 


BGHL89 


Smll 


4577.688 


0.248 


-0.775 


-300 


300 


BGHL89 


Smll 


•4642.232 


0.378 


-0.520 


-150 


150 


BGHL89 


Lu 11 


3819.670 


0.000 


0.510 


-150 


300 


LWDS0I 


En II 


3907.110 


0.207 


0.170 


-150 


200 


LWDS01 


Eu II 


4129.720 


0.000 


0.220 


-150 


400 


LWDS0I 


Eu II 


•4205.040 


0.000 


0.210 


-50 


200 


LWDS01 



References: B BPL89 = Blackwell et al. 119891 . B G80 = Biemont & Godefroid ll980l , BGHL 81 = 
Biemont et al. 179871 . CSSTW82 = Cardon et al. 179821 . DLSC03 = Den Hartog et al. 120011 . G67 
= Gallagher 119671 . GS94 = Gratton & Sneden 119941. HL GBW82 = Hannaford et al. J1982J, HS80 = 
Huber & Sandema n [79801 . LBS01 = Lawler et al. 120011 . LWDS01 = Lawler et a l. tXKm . MB87 = 
Migdalek & Baylis 119871, NKWL99 = Nitz et al. 119991 . PBL95 = Pinnington et al. 119951 . VALD95 = 
see Piskunov et al. 119951 



Table B.l. Assigned average values of eri og g f for each element. 



Element 


^log 9f 


Element 






0.10 


Ni 


0.03 


Ms 


0.07 


Zn 


0.10 


Al 


0.11 


Sr 


0.10 


Ca 


0.11 


Y 


0.03 


Sc 


0.04 


Zr 


0.03 


Ti 


0.05 


Ba 


0.03 


V 


0.05 


La 


0.03 


Cr 


0.05 


Cc 


0.10 


Mn 


0.06 


Nd 


0.03 


Fe 


0.03 


Sm 


0.05 


Co 


0.10 


Eu 


0.03 



